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We apply a general method developed recently for the derivation of the diagonal 
representation of an arbitrary matrix valued quantum Hamiltonian to the particular 
case of Bloch electrons in an external electromagnetic field. We find the diagonal 
representation as a series expansion to the second order in h. This result is the 
basis for the determination of the effective in-band Hamiltonian of interacting Bloch 
electrons living in different energy bands. Indeed, the description of effects such as 
magnetic moment-moment interactions mediated by the magnetic part of the full 
electromagnetic interaction requires a computation to second order in h. It is found 
that the electronic current is made of two contributions: the first one comes from the 
velocity and the second one is a magnetic moment current similar to the spin current 
for Dirac particles. This last contribution is responsible for the interaction between 
magnetic moments similarly to the spin-spin interaction in the Breit Hamiltonian 
for Dirac electrons in interaction. 

PACS numbers: 



I. INTRODUCTION: 

The properties of electrons in solids are usually described in the framework of Bloch theory 
of electrons in a periodic potential. In particular, the study of the dynamics and transport 
properties of Bloch electrons perturbated by external fields led to important results for the 
understanding of metals, semiconductors and insulators properties |l| . For weak fields such 
that interband transitions called Zener tunnelling are negligible, the dynamics of a Bloch 
electron in a given nth band is usually based on the following semiclassical equations of 
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motion 

r = d£(k)/hdk 
hk = -eE - er x B(r) (1) 

where E and B are the electric and magnetic fields respectively and S(k) = £o(k) — m(k).B 
is the energy of the band including a correction due to the orbital magnetic moment m(k) . 
In the band energy £q of the unperturbated crystal, the electron momentum K has been 
substituted by the gauge covariant momentum k = K + eA(R)/h. This substitution has 
been first justified by Peierls in the context of the tight binding model and for this reason is 



called the Peierls substitution [2]. A full justification was later given by Kohn [3|. A simpler 
version of the proof was later provided by Blount, Roth and Wannier and Fredkin who could 
derived an approximate band energy operator as an asymptotic series expansion in the fields 
strength as well as general expressions for the first few terms in this series {4] {5] 6]. Their 
methods, although different, are all based on some approximate unitary transformation 
of the initial Hamiltonian of Bloch electrons in an electromagnetic field which eliminates 
the interband matrix elements and leads to an effective diagonal-in-band energy operator. 
The principal advantage of the Blount's method is that it is an application of a general 
diagonalization scheme applicable to any kind of matrix valued Hamiltonian whereas the 
two other ones are specific to the solid states. In particular, Blount has also considered 
the case of a Dirac particle in an electromagnetic field whose Hamiltonian was diagonalized 
to the second order in the fields strength and first order in their first derivatives. Later 
on, Weigert and Littlejohn developer a systematic method to diagonalize general quantum 
Hamiltonian in a series expansion in h [3] instead in fields strength, an obvious advantage 
of this kind of expansion is obviously that it can be valid for strong external field but 
also that the semiclassical limit is readily obtained. Unfortunately, the method [7] involves 
formal series expansion in terms of symbols of operators which makes the method very 
complicated for practical applications. It is worth mentioning that recently a variant of the 
Foldy Wouthuysen transformation valid for strong fields and based also on an expansion in 
h of the Dirac Hamiltonian was presented |8| . 

This is not the end story with regards to the Bloch equations (OQ). Indeed, Karplus, Lut- 
tinger and Kohn [9] predicted very early a spontaneous Hall effect in ferromagnetic materials 
due to a corrective term to the velocity in Eq. ([!]), known as the anomalous velocity. Later, 



Adams and Blount [10J, by interpreting this term as resulting from the noncommutativity of 
the intraband coordinate operators, derived new semiclassical equations of motion for Bloch 
electrons with an anomalous velocity. However these equations turn out to be correct only 
for external electric fields. It is only recently that the correct equations of motion of Bloch 
electrons in the presence of both electric and magnetic fields and including anomalous veloc- 



ity were derived by Chang and Niu [llj. Indeed, using a time-dependent variational principle 
in a Lagrangian formulation and a description of the electron in terms of wave packets, they 
found the following new equations of motion in the presence of electromagnetic fields 

r = dS{k)/hdk-kx 9(k) 
hk = -eE - er x B(r) (2) 

The correction term to the velocity — k x 6 is the anomalous velocity which is due to 
the presence of a Berry curvature 0(k) of electronic Bloch state in the given nth band, 
associated to the electron motion in the nth energy band. For crystals with simultaneous 
time-reversal and spatial inversion symmetry, the Berry curvature and the magnetic moment 
vanish identically throughout the Brillouin zone. This is the case of most applications in 
solid state physics, but there are situations where these symmetries are not simultaneously 
present like in GaAs where inversion symmetry is broken or in ferromagnets which break 
time reversal symmetries. In the same way, the presence of a strong magnetic field, the 
magnetic bloch bands corresponding to the unperturbated system breaks the time inversion 
symmetries. In all these cases, the dynamical and transport properties must be described by 
the full equations of motion given by Eq. (T2]). Notice that even for crystals with simultaneous 
time-reversal and spatial inversion symmetry, energy bands degeneracies can lead to a non 
vanishing curvature and magnetization; a typical example is provided by the graphene [12I ] 
This case is due to the presence of a topological Berry phase associated to bands degeneracies 
[13I ] . This particular situation is not considered in this paper. 



The method developed in [H] is by construction limited to the semiclassical level, but 
the description of phenomena such as the electromagnetic interaction of Bloch electrons, 
as will be discussed in this paper, requires a theory which goes beyond the semiclassical 
approximation. Recently we came back to the initial considerations of Blount and others 
with regards to the diagonalization procedure for an arbitrary matrix valued Hamiltonian 
H (applicable to any kind of quantum system which has an energy band spectrum) in the 
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presence of external fields and a first original method based on a differential equation of 
the diagonal in-band energy operator with respect to Planck constant h was proposed in 
[jjj]. In this approach, where h is promoted as a (formal) running parameter noted a, 
one has to diagonalize H at the scale a where it is assumed that the canonical dynamical 
operators satisfy the algebra [R l a ,P^\ = iaSij. Relating two diagonalization processes for 
close values of a, leads to a differential equation of the required diagonal Hamiltonian e a with 
respect to a. This differential equation has to be supplemented by an additional equation 
which is the consequence of the unitarity condition of the matrix U a diagonalizing H at the 
scale a. The resolution of this differential equation can then be performed by a systematic 
series expansion in h, and in this way, at least in principle, an exact diagonalization of 
arbitrary Hamiltonians can be achieved. This approach reveals that the diagonal energy 
operator is most naturally expressible in terms of covariant (noncanonical ) coordinates 
r = R+j4r and momentum operators p = P+ylp which are both corrected by Berry 
connections terms ^4r/p and which satisfy a non- commutative algebra. Particle motion in 
this noncanonical phase space is obviously drastically modified (as in Eq. (J2J)). Particularly 
interesting is the fact that, in the semiclassical limit, which is often enough to get physical 
insight to the problem considered, the diagonal Hamiltonian is obtained by a straightforward 
integration of this differential equation. In this limiting case, the differential approach turns 
out to be so powerful that the general diagonal representation for an arbitrary matrix valued 
Hamiltonian in terms of covariant operators and commutators between Berry connections 
could be givenfactually, this general formula was first derived by a direct diagonalization 
procedure in 15|] ) . This result allowed us to deduce effective semiclassical Hamiltonians and 
to predict new phenomena in various physical situations. First, the study of Bloch electron 
in magnetic Bloch bands showed that besides the position operator which get a Berry- 
phase contribution (as already shown by 



111]), the momentum in the band energy £ (k) also 
has to be replaced by a new Berry- dependent momentum operator k = K + eA jfi — A? 
instead of the Peierls substitution. It turns out that this result is essential for the correct 
derivation of the full equations of motion for Bloch electrons Eq. (T2]). Likewise, for electrons 
in graphene in a magnetic field, it was also observed that it is in terms of k that the 
semiclassical quantification of the orbit has to be achieved 171 ]. 

At the semiclassical level other systems were also investigated, like Dirac electrons in 
electromagnetic [14] and gravitational fields with the discovery of a spin-magnetotorsion 
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coupling. The study of the photon in a static gravitational field when polarization effect are 
taken into account predicts the gravitational birefringence phenomenon where an helicity 
dependent anomalous velocity deviates the photon from the usual Einsteinien geodesies 19]. 
This kind of polarization effects which are called spin Hall effect of light have been recently 
observed 20J]- Despite these results, the applicability of the differential approach beyond 
the semiclassical turns out to be very complicated. Even the deduction of the diagonal 
representation of a generic Hamiltonian at second order in h is a prohibitively difficult 
problem, although a solution was found for simple practical applications like a photon in an 
homogeneous isotropic media and Bloch electron in uniform electric field jljj]. But clearly, 
one can not expect to use this method in the case of several Bloch electrons in electromagnetic 
interaction, a problem we would like to consider in this paper. 

However, very recently, a new general and powerful method for the diagonalization of 
an arbitrary matrix valued Hamiltonian has been proposed by Gosselin and Mohrbach (re- 
ferred as GM) 2l|.It leads to a particularly compact and elegant exact expression for the 
required diagonal energy operator. This approach is therefore particularly well adapted to 
problems in solid state physics and to Dirac particles in external fields. This last case can be 
considered intuitively as a simple two bands versions (particles and anti-particles) of Bloch 
electrons in a crystal. The philosophy behind this approach consists in mapping the initial 
quantum system to a classical one which can be diagonalized and then to return to the 
full quantum system. This method is not based on a differential equation for the diagonal 
energy operator with respect to h, but it also requires the introduction of new mathematical 
objects like non-commuting operators which evolve with h promoted as a running variable. 
This new mathematical construction leads us to define a differential calculus on a non- 
commutative space showing some similarities with the stochastic calculus as both stress the 
role of second order terms. This approach allows us to write both the diagonal Hamilto- 
nian e (x) and the transforming matrix U (x) (where x = (r, p) is the phase space of the 
covariant dynamical operators) as a result of the application of integro-differential opera- 
tors on e$ (X ) and Uq (Xq) respectively, i.e. e (x) = O (eq (Xo)) and U (x) = iV (Uq (X )). 
Matrices with the subscribe correspond to the operators replaced by classical commuting 
variables X = (Ro, Pq) • The only requirement of the method is the knowledge of Uq (X ) 
at a = which gives the diagonal form e (^-o) • Generally, these equations do not allow 
to find directly e (x), U (x), however, they allow us to produce the solutions for e (x) and 
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U (x) recursively in a series expansion in H. But contrary to the procedure in 0], it turns 
out that the expansion in a series of fi is much more easier to obtain than by the succes- 
sive integration of the differential equation and therefore more convenient for getting higher 
order contributions. Remarkably, it was also found that the exact expression for e (x) is 
actually an exact solution of the differential equation of [14]. This result obviously places 



the present approach on a firm base|2l|]. Another interesting feature of this approach is the 
confirmation of the fundamental role played by Berry curvatures in these systems since the 
method results in an effective diagonal Hamiltonian with Berry phase corrections as well as 
noncommutative coordinates and momentum covariant operators as in previous approaches 
QQ). 

Although similar in spirit to Blount's method and in particular to Weigert and Littlejohn 
one Q, the approach proposed in 0] is essentially different as it is based on a very new 
mathematical formulation. (The general method of p| leads also to a diagonal in-band 
energy representation as a formal series expansion written in terms of symbols of operators 
which makes the method very complicated for practical applications). In our opinion this 
new approach is more tractable for applications and, as an illustration of this statement, a 
general in-band energy for any arbitrary Hamiltonian to the second order in h was achieved 
in 2l) (higher order expressions becomes again very cumbersome but could in principle be 
computed). This expression will be the starting point for a straightforward study of a single 
Bloch electron in an external electromagnetic field. Another purpose of the present work 
is the adaptation of GM's results to the case of Bloch electrons in interaction. Note that 
both problems can be transposed to the case of Dirac electrons which are actually treated 



by the same method in the another paper 22j. It is obviously the Coulomb (electrostatic) 



interaction between Bloch (Dirac) electrons that dominates over the magnetic one, so that 
a first order diagonalization seems to be sufficient. But, in the presence of non vanishing 
electronic magnetic moments, other effects like moment-moment interactions mediated by 
the magnetic part of the full electromagnetic self-interaction are expected. This comes out 
by analogy with the spin-spin interaction in the Breit Hamiltonian of non-relativistic Dirac 



particles 23| which is recover in the non-relativistic limit [22]. Obviously these kind of 
interactions are of second order in h and a diagonalization procedure which goes to this 
order is necessary. 

Before starting, two points are worth mentioning. As already noted in Q] but also in 
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[ijj] [ijj], there is certain latitude in finding U which reflects a kind of gauge invariance of 
the method. Because of this freedom two equivalent diagonalizing operators leads to two 
different forms for the diagonal Hamiltonian but of course to the same eigenvalues. This is 
similar to the Schrodinger equation in a magnetic field where the Hamiltonian and the wave 
function gauge dependence combines itself to give gauge independent energy levels. In our 
case it turns out that this freedom is only present at the second order in h, because the first 
order diagonalization is performed with the zero order matrix Uq (X ) which can be uniquely 
defined. Actually this gauge dependence can be included in the gauge covariant dynamical 
operator x = (r,p). Then the diagonal in-bands energy operator is uniquely defined when 
it is written in terms of x instead of the canonical operators X. As a particular gauge choice 

nn 

can be made on the ground of simplicity and convenience, here as in previous works 14 21 
the reality condition of diagonal elements (the anti-hermitian diagonal elements are setting 
to zero) of U is imposed. 

The second point we would like to mention is that, in order to simplify the expressions, 
only time independent electromagnetic fields are considered, but results can be easily ex- 
tended to include time dependent interactions. In addition, as the ultimate goal is to consider 
Bloch electrons in interaction through an internal electromagnetic field, the time dependence 
of the vector potential can be safely neglected as all retarded effects in the electromagnetic 
interaction are negligible owing to the fact that Bloch electrons are non-relativistic. 

The paper is organized as follows. In the next section we give a reminder of the diago- 
nalization procedure of [^j] for an arbitrary matrix valued Hamiltonian. We provide some 
detailed formulas for the diagonalized energy operator at the second order in h as well as 
for the Berry phases at this order. Section 3 applies this formalism to the case of a Bloch 
electron in an external electromagnetic field. The diagonalized Hamiltonian at the second 
order in H is written in terms of the transformed dynamical variables and magnetization 
operators. In section 4, we consider the case of P Bloch electrons interacting through an 
internal electromagnetic field. Diagonalizing the matter part of the Hamiltonian to the 
second order in h and solving for the electromagnetic field yields the P particles effective 
Hamiltonian. Last section is for the conclusion. 
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II. DIAGONALIZATION OF AN ARBITRARY MATRIX VALUED 

HAMILTON! AN 

To start with, an outline of the approach developed by Gosselin and Mohrbach to diago- 
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ize formally an arbitrary matrix valued Hamiltonian is given, with the notations of GM 
2l|. Let consider an arbitrary quantum mechanical system whose state space is a tensor 
product L 2 {K Z )®V with V some internal space. In other words, the Hamiltonian of this 
system can be written as a matrix H (R, P) of size dimV whose elements are operators 
depending on a couple of canonical variables X = (R, P) . The archetype example is usually 



the Dirac Hamiltonian with V = C 4 , but as shown in 15( the following set up fits with all 
system presenting an energy band spectrum, as for Bloch electron (where V correspond to 
the energy band indices), a system which is the main concern of the present paper. 



In 2l[, a method to find an unitary matrix U (X) to diagonalize any arbitrary matrix- 
valued quantum Hamiltonian H (X) such that e (X) = UH (X) U + is the diagonal in-band 
energy operator was achieved. 

As explained in the introduction, the principle of this method is to link continuously 
an Hamiltonian in which the variables are considered as classical (i.e. h = 0) to the true 
Hamiltonian we aim at diagonalizing (that is at scale K). The idea is to diagonalize the 
Hamiltonian for ft = 0, which appears in general to be much easier, and then to come back 
to scale h to obtain the required Band Hamiltonian. To perform this program, we have to 
proceed in an indirect way. Actually, we first need to introduce a family of canonical variables 
(R a , P Q ) indexed by a continuous parameter a G [0, h], such that the commutators are given 
by [R^, Py = ia and then, to introduce for any arbitrary function F (R a , P Q ) both a notion 
of differentiation and integration describing the variation of F (R Q , P a ) as a varies. These 
notions have to take into account the fact that the commutation relations are depending on 
a. It leads us naturally to introduce the notion of infinitesimal non commutative canonical 
variables (dH a , dP a ) as well as a form of differential calculus presenting some formal analogy 
with the non commutative stochastic calculus. The introduction of these differentials notion, 
will allow to connect ultimately our Hamiltonian at scales and h. 

More precisely now, we introduce a space of non commuting infinitesimal operators 
dX l a = {dR} a1 dP^} i = 1,2,3 indexed by a continuous parameter a, that satisfy the fol- 
lowing infinitesimal Heisenberg algebra with a reversed sign [dR l a , dP^,,] = —ida5 a , a ,5ij and 
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[dR l a , dR 3 a ^\ = \dPa, dPa,] = 0. From it, we define a set of running coordinate and momen- 
tum operators by writing the following formal sums R l a = R 1 — dR\, and P l a = P % — J h dP\ 
with the choice of convention dR l a = R l a — R l a _ da and dP % a = P % a — P % a _ dOL , so that the running 
operators satisfy [R l a , P^] = ia5ij and [R l a ,R J a ] = [P^,P^] = 0. For a = h we recover the 
usual canonical operators R l = i?^and P % = P % h which evidently satisfy the canonical Heisen- 
berg algebra. The differential of an arbitrary function F(X. a ,a) where X l a = {i? Q ,P^} on 
this space is given by 



6 ,6 

dF (X a , a) = V ^ (X a , a) dX\ - - V x * V x , F (X Q , a) (dX^dX^ + dX l a dXi) 

+ +(F(X a> a)))«fa (3) 

with i,j = 1..6. We also assume that X l a = R % a for i = 1,2,3 and X l a = P % a for 
i = 4,5,6 .The notation (F (X a ,a)) (which in [14j was corresponding to the operation 
— lAsymV^ VpiF (X a , a)) is defined as a specific procedure on a series expansion of F 
in the variables R % al P % a in the following way : let F be a sum of monomials of the kind 
Mi (Rq,) M 2 (P q ) M 3 (R a ) .... the Mj being arbitrary monomials in R a or P a alternatively. 
Let the operator V^Vpi acts on such an expression by deriving all combinations of one 
monomial in R a and one monomial in P a . For each of these combinations, insert a dR l a at 
the place where the derivative Vj^ is acting and in a same manner a dP^ at the place where 
the derivative Vpi is acting. This leads to an expression with two kind of terms, one kind 
being proportional to the dR l a dP^ and the second proportional to dP^dR l a . Then rewrite 
this expression in terms of dR l a dP^ + dP^dR^ and dR l a dP{ — dP^dR l a = —i5^da. Then 
(F (X a ,a)) is defined as minus the contributions of terms proportional to —iS^da in the 
computation in the procedure just considered. This definition implies a procedure which is 
clearly dependent of the symmetrization chosen for the expansion of F. 

To make the definition of (F(X Q ,a)) clearer, consider some important practical exam- 
ples. If the function F has the following form F = \ {A (R Q ) B (P Q ) + B (P Q ) A (R a )) 
which corresponds to a frequent choice of symmetrization in R a and P a , then (F (X a , a)) = 
| [A (R Q ) , B (P a )]. Another choice of symmetrization leads in general to a different result. 
For instance, if we rewrite the same function F in a fully symmetrized form in R Q and P Q 
(that is invariant by all permutations in R a and P a ) which is also often used, we have now 
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have a different result since (F (X a , a)) = 0. 

Nevertheless, this dependence of (F(X. a ,ct)) in the symmetrization choice is not aston- 
ishing at all. Actually changing the symmetrization of a function F (X Q , a) introduces 
some explicit terms in a which changes also the term d a Fda present in the differential Eq. 
(J3J). As a consequence, neither the partial derivative with respect to a, nor the bracket are 
invariant by a change of form. But, what is invariant is the sum d a F + (F). This assertion 
is shown in 21]. 

Now, what really matters for us in Eq. ([3]) is this invariant term proportional to da this 
is why we define an expectation operation £ (.) so that 

£ (dF (Xa, a))=s( ( 9F( f^ a) + (F (X a , a))) da) (4) 

Combining the expectation operator as well as the differential allows to relate a function 
evaluated at the physical scale H, F (X^, h) to that same function evaluated at the scale 0, 
F (Xo, 0) that is when the canonical variables are considered as classical ones. Actually by 
integration of the previous relation one has: 

F (X,, h) = £^F (X , 0) + J ( ( ^^^ + ( F ( X - «))) da ) } (5) 

(we use a property that the expectation operator satisfies £ (F (X&, h)) — F (X&, h)). 

Here recall that we denote by X Q the dynamical variables when P and R are considered 
as commuting classical variables that is when h = a. This last formula is the starting point 
of the method. Iterating this relation, by successive differentiations and integrations, one 
can relate a full quantum function F (X^, H) to the same function evaluated with classical 
variables X . This will prove very useful in a problem of diagonalization of a matricial 
Hamiltonian since in that case, the diagonalization when the canonical variables Xo commute 
reduces to the diagonalization of an usual matrix of finite size. Thus diagonalizing our 
Hamiltonian at scale 0, that is finding an F (X , 0) is in general an easier task. 

However, this not the end of the story. Having found a way to relate F (X^, K) to its 
classical counterpart is not enough since we want ultimately to recover expressions of interest 
evaluated at X^ not at the fi = scale. This kind of coming back process after a " classical" 
diagonalization is performed by an other operation defined in {2]]]. It allows in the previous 
integral relation to replace, inside the expectation, X and X Q by X^ at the price of a 
modification of the expression inside the integral. Define the exponentiated Bracket plus 



11 



Shift operator (EBS) between a 2 and a.\ as acting on any function F(X Q2 ,a 2 ) to yield an 
other function depending on (X ai ,a 2 ), a i > a 2 '■ 



exp (-(■)f 2 _ >Ql ) = Texp 



Sx ai (-) a S Xa da 



a-2 



s 



Sy 



(•)/?! ft 



/3l 



d(5i...d(3 n 



'a2</3n<.../3i<ai 

where the Shift operation Sx ai sets the dynamical variables X Q to X ai and satisfies 



iS*x„ Sy 



SyL a whatever the values of a and f3. Apart from the repeated application 



of the Bracket (.) the EBS operation has the virtue to shift progressively the variables 
from X , the "classical variables" to X = X^ the full quantum variables. Actually its main 
property is the following : 

£F (X Q2 , a 2 ) = £ exp (- (.)f 2 _ >Q1 ) F (X aa , a 2 ) 

so that it can of course be specialized to : 

SF (X , 0)=E exp (- (.> ^) F (X , 0) 

These two formulas can be understood intuitively as follows. The EBS operation changes 
the function (by the action of the bracket defined above) but also changes progressively 
the variables from X a2 to X m (through the shift operator). As a consequence, and despite 
the appearances, the expression in the right hand side £exp ^— (-)f 2 ^ Ql j F (X Q2 , a 2 ) is a 
function of X Ql . Moreover, both the EBS operation and the shift of variable compensate 
each over to produce the equality with the left hand side. 

With this mathematical construction in hand it was possible to write the solution of 
our diagonalization procedure for a general matrix valued Hamiltonian through an unitary 
transformation U as the solution of the following system of integro-differential equations. 
Introducing Uq (X ) the diagonalization matrix when h — 0, such that the classical energy 
matrix obtained as Eq (X ) = UqH (X ) Uq is a diagonal matrix, we could write : 



e(X) 
U(X) 



£ 



T exp 
Texp 



0<a<h 



0<a<h 



£ o (Xo 
^o(Xo) 



(6) 
(7) 



T is the notation for the time ordered product, the operator exp ^— (•)f 2 _ >Q , 1 j acts as ex- 
plained above, and O a and N a act in the following way : 



12 



The operation O a e (X a ) = + (.)) e a (X Q ) can also be written as O a eo (X a ) = 
(T a + M a ) e (X a ) with a translation operator T and a "magnetization" M operator (this 
terminology is explained in [ijj]) where 



T a e (X a ) = |^ (A%V Rl e (X a ) + V^o (X Q ) ^ + ^V Pi e (X a ) + V Pl e (X a ) A 
M a e (X a ) = [e (X a ) , A*] A% - [e (X a ) , ^] ^ } + H.C.} 



+V 4 



These equations require some explanations. First, the operators A^ are given by 



iP 



and A a (X a ) 



[C/ Q (X a )VpC/+(X Q )] 
-i [C/ a (X a )V R [/+(X Q 



(9) 



and V+ and "P_ are respectively the projectors on the diagonal and off-diagonal elements of 
matrices. From Eq. ([8]) one can deduce the following relation 



8 ( Texp 



T n da 



0<a<h 



£o (Xc 



£e (x) 



where x = (r, p) and r and p are new covariant coordinate and momentum operators defined 
in the following way : 



r = R+,4 R 



P 



(10) 



The Berry connections terms being defined as : 



V + [AS] da+ [ \\\v + [A?] • V x / V + [A*] 

t<h JQ<a<h Z LL J0<a 1 <a 

A P = ! V + [Al] da + I \\\v + [Al] . V X / V + [<] 

J0<a<h J0<a<h^ LL J0<ai<a 



+ H.C. 
+ H.C. 



da\da + 
daida + 



and X denotes the vector (R, P). In fact in 21] we show that A n and A p have to be 
corrected by a third order in h terms 5 A 11 and 5A F but these terms will always be neglected 
here. 

Therefore, the operator T as a part of O naturally leads to the emergence of the covariant 

re be the physical dynamical variables of particles 



dynamical coordinates which turn out t 
as shown in several situations [ijljljj] 19] [3] 
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Second, using the gauge setting to zero the anti-hermitian diagonal elements of U a [21] 
we have 



N a U a 



+ l -V-{[E a ,A*]J*- [s a ,A5]A*}+H.c] - % -V-{[A*,J*]U a } 



where the inverse of the commutator operation [, e a ] has the following properties 

[[.,e a ] _1 .M,e a ] = [.,e a ] _1 .[M > e a ] = M for [M,e a ]^0 

Ue a \-\M = if [Af,e a ] = (12) 

for an arbitrary matrix valued operator M. It means that the operator [.,£ a ] _1 acts on the 
space of endomorphism of the state space in the following way : it is zero when acting on 
the kernel of the operator [., e a ] whose action is to compute the commutator with e a , and is 
the inverse of [.,£«] on the complementary subspace of the kernel. 

Note that, as it will appear clearly later on, having both e (X) and U (X) at order n in 
h, and reinserting in the exponential of Eqs. ([6]) ([7]) allows us to find e (X) and U (X) at 
order n + 1 in H. A needed assumption for this procedure to work is that the diagonalization 
Sq (Xq) = UqHq (Xq) Uq is explicitly known when h = 0, i.e. when P and R are treated as 
commuting variables. 

We end up this section with a technical remark that will be important for the sequel. As 
shown in our formula (JTj) the solutions for the diagonalization process depends on e (Xo), 
Uq (X ) . While the final results do not depend on the way variables are symmetrized (that 
is the order we write the products of components of X, one has to start with an initial 
symmetrization for the diagonalized energy at the zeroth order as well as for Uq (X ) or 
equivalently the Berry phase. Since this detail will be important only at the second order in 
h while considering the Bloch electron in an electromagnetic field, we do not mention any 
choice for the moment. 

Eqs. (J7J) (jlip show that a diagonal Hamiltonian representation can be found to any 
desired order in h, and in the following we will carry it out until the second order. Let start 
first with the first order. 
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A. Covariant dynamical operator algebra 



From equations Eq. ({TO]) we readily deduce the following non trivial algebra between the 
dynamical operators 

[r u rj] = ih 2 e% = th 2 (V Pl A R] - V P] A Ri ) + h 2 [A Rj ,A Ri ] 
\p i7Pj ] = ih 2 ®% = -ih 2 (V Rl A Pj - V R .A Pi ) + h 2 [A Pi , A Pj ] 



-thdij - ih 2 {V R A Rj + V Pj A Pi ) + h 2 [A Pi ,A Rj ] (13) 



where the terms Gjj are Berry curvatures definitions. Of course these non trivial commuta- 
tion relations also give new contributions to the equations of motion and thus lead to new 
phenomena Jjj] 19] 18]. The commutation relations are valid to any order in h, but 

in practice we can compute them as well as the energy e (X) in a series expansion in h. 



B. Diagonal representation at the first order in h 

In this section we derive by a straightforward application of Eq. ([6]) the semiclassical 
effective diagonal in-bands Hamiltonian e for one particle. Using the fact that at this level 
Jq O a dae (X a ) = h (T h + M h ) e (X) we have 

e = so (X) + h (T h + M n ) e (X) + 0(h 2 ) (14) 

From the definition of T and M in Eq. ([8]) we see that we must determine the quantities 
^4^" = (A^,Af t ) Eq. ([9]) defined previously. At this level of the approximation and due to 
the factor h in the previous equation it is enough to know this quantity at the zeroth order 
in H. Skipping the h index, they are thus simply given by 



Af = -4 R = zf/o(X)V P f/ + (X) 
Al = A* = -iU (X) V R f/+ (X) 

where as before Uq is the matrix that would diagonalize the Hamiltonian if R and P where 
commuting variables (that is the matrix Uq which diagonalizes the classical energy when 
h = 0). It means that, interestingly, at this first order in H, the diagonalizing matrix U is 
the matrix defined at the order h° in which the classical variables X have been replaced by 
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the quantum operator X. This is the clue of the method: at each order H n we just need to 
know the matrix U at the order If 1 ' 1 . Therefore, the semiclassical diagonalization is within 
the present approach simple to achieve (and simpler than any other method we are aware 
for) as it only requires the determination of the matrix diagonalizing the classical energy. 

Introducing also the covariant dynamical operators x = X+v4 x with the Berry con- 
nection j4 x = hV+ [A^~\ , as well as the integrated non projected Berry connections 
A x = Jq A^da = ^L4 X at this order, the diagonal energy Eq. (TTJJ can be written as 



e = eo (x) + { [e„ (X) , A*] A p > - [e Q (X) , A*] A* } + 0(h 2 



(15) 



where summation over i is assumed. Here eo (x) corresponds to the classical diagonal energy 
operators in which the classical variables Xo have been replaced by the full noncomutative 
covariant operators which satisfy a non commutative algebra as we will see later when 
considering the practical examples. The second contribution which is of order H, gives for 
instance for Bloch electrons in a magnetic field B as we shall see later, the coupling between 
B and the magnetic moment operator . 

Note also that this general diagonal energy for an arbitrary one particle system was 
already derived in previous works by different methods 15] 14 and turns out to be very 
useful for the study of several different physical systems lfjf 18( 19( 171 ]. 



C. Diagonal representation at the second order in H 

The Hamiltonian diagonalization at this order requires the matrix U at the first order: 
U (X) = U (X) + HUx (X) U (X) where U x (X) is determined from the relation HU X (X) = 
Jq N a daUo (X) as a consequence of Eq. ([Zj) . With expression (fTTj) we readily obtain : 



^(X) = [.^o]" 1 



V -{1 C^o^o ( X ) + ( X ) ^ + A o V ^ ( X ) + V ^o (X) A$) } 



e (X),A* 1 



A Pl - 



*/-\.r. 



A Rl A Pl 



(16) 



At the same order the (non-diagonal) Berry connections ^4 X = (A^, A^ ) are again given by 



A* (X Q ) = i [U a (X a ) V P C/+ (Xa)] = -U a (X„) R a U+ (X„) 

a 

and Al (X Q ) = -i [U a (X Q ) V R U+ (X q )] = -U a (X a ) P a U+ (X a ) 

a 



(17) 
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where now U (X a ) is the transformation to the first order in a i.e. Uq (X a ) + aU\ (X a ) , in 
which X is replaced by the running operator X a . Using the hermiticity of A*, so that one 
has A* = A* + (A£) + we can expand A x as : 

A* = Q ([1 + aU x (Xa)] U (Xa)) ^ (ttf (X a ) [1 + aUt (X a )]) + ijET.<7.^ - ^ 

(the - factor reminds that in our definition of .4.* the Gradient with respect to X a is 

normalized, i.e. divided by a). After some recombinations, the previous expression can be 

written in a more convenient form : 
1 



A 



U (X a ) [X a ,U+(X a )]+H.C. 



2a 

+ l - [U x (X Q ) [U (X a ) X a U+ (X a )] + [U (X a ) X Q [/ + (X a )] [/+ (X a )] + ff.C. - X c 

using now the fact that at the lowest order Uq (X a ) [—, Uq (X a )l = Aq (X a ), one has : 
1 



At 



2a 
1 



[U (X a ) [X a ,[/ + (Xa)] + ff.C.] 



+ - [C/x (Xa) Af (Xa) + ^ (Xa) ^ (X a ) + [X Q , ^ + (X Q )] + H.C.] 



decomposing (X a ) into Hermitian and antihermitian part we are thus led to : 
A* = (£/ (X Q ) [X Q ,f/+ (X a )] + ff.C.) + [X a +^,a/i (C/+ (X a ))] 

+ (AfH (U+ (Xa)) + ft (f/+ (X Q )) .4*) 

where ah (Z) and a/i (Z) denote the antihermitian and Hermitian part of an operator Z 
respectively. Now with Eq. ffl6|) and after the integration, we are led for ^4 X = HA$~ + h 2 A* 
to the following expression : 

= *Z (R+f P+f X 1 ) - [B, X Q + Af] 

where we introduced the notations 



(18) 



A* (r+| J®, P+f A?) =Af + \ {^VnX + < l ^nX + B.C. 



(19) 



and 



B 



1 ■ (v- {\a*> V Rl e (X) + l -A$V Pl e Q (X) + B.C. 
-- { [e (X) , X'] A* - [e (X) , < ! + H.C)\ 
•' £ ° rl • ( p - {1 A * Vr * £o W + ^o P; Vp ; eo (X) + H.cX 



-- <! V-Ao l V + Ao' - V-Ao l V+Ao> + H.C. 



(20) 
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These formula, although abstruse will be useful in our next sections. Having found the 
matrices U (X) and A^ at the required order we can now determine the series expansion for 
the diagonalized Hamiltonian to the second order in h. From equation Eq. ([6]) we can write 



e (X) = e (X) + / O a dae pt a ) + / O ai / O a2 da 2 da 1 e (X a2 ) - - (e (X)) 
Jo Jo Jo ^ 

where the last contribution is given by 

~h (e„ (X)) = % -hAsym {V Pl V Rl e (X)} 
The first contribution J O a dae (X) can be expanded as before as: 

J O a dae {X) = J P+j^^V^eolX^ + V^eolX^^+^Vp^olX^ + Vp^olX,)^^}^ 

+ / V { ~ { h (Xa) , ^] 4? - [*0 (x a ) , ^] ^ } + #.c. 



+ 



da 



where the Berry connections have to be expanded to the first order, whereas the second 
order contribution 



\ O ai / O a2 da 2 da 1 e (X) - - (e (X)) 
Jo Jo * 



has to be expanded to the zeroth order in the Berry connections. Notice that due to the 
integration process, the squared terms in A^, Aq as well as the first order terms in the 
Berry phase get a | factor. 

The consequence is that these contributions can be recombined to yield : 



e(X) = e (x) + -V + \ e (x),A Rl 



A 1 



e (x),i J 



A 



Ri 



So W 



]} 



--V 4 



Bo (X),i^ 



A p '- 



E (x),^ 



A Rl A Rl 



£o (x) , A 



Rt 



with 



A Rl 
A Pl 



A 1 



e (x) , A 1 



A R \A P > 



A 



Rt 



A Pl 
h 



{so (x)> 



(21) 



1 - ~ (V+J$.Vy) 

1 - ~ (M X -Vx) 



A Rl (x) + E.G. 
A Pl (x) + E.G. 
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and where we have denoted again ^4 X (X) = f A* (X) da. Remark that the last term 
— | (eq (x)) in the expression for e (X) is of order h 2 since (sq (x)) involves some commutators 
and is thus of order H. As explained above, one can, at each order of the expansion, safely 
replace the canonical operators X = (R, P) by the covariant ones x = (r, p) which are given 
by the following expression to the second order (dropping once again the h index) : 



h 2 

x = X+v4 x = X+/L4 X + - Af 
with : ^4 X = (^4 R , ^4 P ) and (as before we drop the index K) X 



X, 



A 



R 



0<a<h 

V 



0<a<h 



J r + [jg\da+ I \ \{V + [A*].V* + H.C) J V + [A*] 



0<»i<a 



(22) 
+ H.C. 



da i da 



A 1 



hA* (R+^, P+~ a A^ - | [B, R + .4?]] + | (V + [Af] .V X P + [A*] + H.C.) 
J V + [Al]da+ J \ 1(V + [A*].V X + H.C.) J 7> + [<] 



0<a<h 

V 



0<a<h 



0<cn<o 



+ H.C. 



da\da + 



KAl (r+Jm* P+^) -y[ B ' P + ^] + T ^ ' VxP+ ^ + ^ 
Remark ultimately, that had we chosen the variables X rather than x to express our Hamil- 
tonian, we would have rather written at our order of approximation : 



A Rl 
A p i 



1 + - pVtf.Vx) 
1 + \ pVtf-Vx) 



A Rl (X) + H.C. 



A Pl (X) + H.C. 



However, as explained before, the transformed variables x fit better to write the Hamiltonian 
since they enter directly in e (x). 

Eq (j21l) is the desired series expansion to the second order in h of the diagonal Hamilto- 
nian. In the next section it will be the angular stone for the computation of the effective 
in-bands Hamiltonian of a Bloch electron in an external electromagnetic field. Once this 
will be achieved the case of several Bloch electrons will be investigated. 



III. BLOCH ELECTRON IN AN ELECTROMAGNETIC FIELD 

To start, an outline of the first order diagonalization for the special case of an electron in 
an crystal lattice perturbated by the presence of an external electromagnetic field considered 
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refs. and Q 



m rets. |lb|J and [lo|J is given. 

The Hamiltonian of an electron in an crystal lattice perturbated by the presence of an 
external electromagnetic field is 

H= { 6 ) +V p (R) + eV(R) 
2m 

(e < 0) where V p (R) the periodic potential, A and V the vector and scalar potential respec- 
tively. Our purpose is this section is to compute the diagonal in-bands energy Hamiltonian 
for this system to the second order in a series expansion in h. This can be done by using 
the general results of the previous section, in particular with Eqs. ffTB^) (T2"Tj) . The major 
difficulty to find the diagonal representation relies on the fact that in presence of an electro- 
magnetic field, the lattice translation operators T do not commute any more (see [ll| and 
references therein). To deal with this problem it is convenient to express the total magnetic 
field as the sum of a constant field Bo and small nonuniform part 5B(R). The Hamiltonian 
can be written H = Ho + eV^(R), with Hq the magnetic contribution (V being the electric 
potential) which reads 



H = — (P-eA(R)-e5A(R)) 2 + K(R) (23) 
2m 



where A(R) and 5A(R) are the vectors potential of the homogeneous and inhomogeneous 
magnetic field, respectively. V p (R) is the periodic potential. The large constant part B 
is chosen such that the magnetic flux through a unit cell is a rational fraction of the flux 
quantum h/e. The advantage of such a decomposition is that for 5A(R) = the magnetic 
translation operators are commuting quantities allowing to exactly diagonalize the Hamil- 
tonian and to treat <5A(R) as a small perturbation. The state space of the Bloch electron is 
spanned by the basis vector \n, k) = |k) ® \n) with n corresponding to a band index and k 
a common eigenvalue of the translation operators. In this representation K \n, k) = k \n, k) 
and the position operator is R =id/dK. which implies the canonical commutation relation 
[Rj, K ■] = iSij. Note that here R and K, rather than R and P will play the role of canonical 
variables. 

The diagonalization of the Hamiltonian in Eq.( f23|) is first derived for 6 A = by di- 
agonalizing simultaneously Hq and the magnetic translation operators T. Start with an 



arbitrary basis of eigenvectors of T. As explained in 15|, in this basis H can be seen as 



a square matrix with operators entries and is diagonalized through a unitary matrix U(K) 
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which depends only on K (since U should leave K invariant, i.e., U~KU + = K), such that 
in-bands energy matrix is e = UHU + = £o(K) + V(UHU + ), where £o(K) is the unpertur- 
bated (5 A = 0) diagonal energy matrix made of the magnetic bands elements £o,n(K) with 
n the band index. 

Now, to add a perturbation 5A(R) that breaks the translational symmetry, we have to 
replace K in all expressions by 

n = K - |<L4(R) (24) 

the (band) electron momentum and as the flux 5B on a plaquette is not a rational multiple 
of the flux quantum, we cannot diagonalize simultaneously its components Ilj since they do 
not commute anymore. Actually 

ftflP, IF] = iee ijk 5B k (K) (25) 

To deal with this non-commutativity, we adapt our method to diagonalize the Hamiltonian 
perturbatively in h. To start, an outline of the first order diagonalization for the special case 
of an electron in an crystal lattice perturbated by the presence of an external electromagnetic 



field considered in refs. [16) and 15) is now given. 



A. Semiclassical diagonalization: Generalized Peierls substitution 

Following section II, the diagonalization at the lowest order is just obtained by replacing 
U (K) by U (II) . This last matrix would actually diagonalize the Hamiltonian if R in 5A(R) 
was a parameter commuting with K). Note that a subtlety arises here (and that we will 
find again later) in the application in our method. We do not consider, at this level that 
R and K commute, but only that the R in <L4(R) and V (R) commutes with K. In other 
words we assume that R has been replaced by a parameter in these potentials. The reason 
of this difference with our general set up comes from the fact that the initial diagonalization 
is not performed for a function of, say, P alone, but both of P and R through the periodic 
potential. However this difference does not alter our method which allows to recover the 
contributions of the electromagnetic potential as a series of H. 

As consequence of our procedure, the non projected Berry connections are Ajf* = 
iUV KiU + and Aq 1 = eV R l 5A k (R)AQ k . However, it turns out be more relevant to replace K 
by the covariant momentum II in the physical expressions so that instead of Aq 1 we better 



21 



consider the quantity A^ 1 = eA$ k B m e krn i. Remark that the Berry connections just defined 
are non diagonal and are matrices whose index correspond to interband transitions. 

The physical dynamical variables for the ra-th band dynamics to the first order in h 
imply a projection on the n-th band. For the intraband coordinate operator r = r n = 
V n (U(U)RU + (U)) we obtain 

r~~R+Af(7r) + 0(H) (26) 

and for covariant intraband momentum we obtain in the same manner ir = ir n = II + Ajf 
which writes also 

K~n + eA*(Tr)x8B(r)/h + 0(h) (27) 

with = HV n (A^) and V n the projection on the n-th band. Remember that previously V+ 
was the projection on the diagonal elements of a matrix. Using now the general expression 
for the semiclassical Hamiltonian Eq. ffl5l) . we obtain the desired semiclassical nth-band 
Hamiltonian (dropping the index n) e (tt, r) as : 

e (tt, r) = e (vr) + eV (r) - M(-K).5B(r) (28) 

where Eq (it) is the unperturbated nth-magnetic band energy in which K has been replaced 
by 7r, a procedure that we can adequately call the generalized Peierls substitution and which 



was introduced for the first time in 



161 ] . The second term in Eq. ( 128]) constitutes the 



electrostatic potential, and the third is the coupling between <5B and the magnetic moment 
(or magnetization) which is given: 

M(n)=V n (— [e (7r),M?(7r)] x M r (tt)) 



It can explicitly check that this expression of t 



expression found with different approaches 



le magnetization is the same than previous 



24]. 



1. Dynamical operators algebra and equations of motion 

From the dynamical operators a new algebra has to be considered . Indeed we have 

[r\ r j ] = iQ{n) ij 
n[7r\7r J ] = iee ijk SB k (r) + ie 2 e ipk 5B k e jql 5B l Q pq /K 
[r\hTr j ] = ihS ij -iee jlk 5B k (v)e u (n) (29) 
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with Q^(tt) = d l A 3 (7r) — &>A % (7r) the Berry curvature. Whereas the term of order 6B 2 in 
[ir l , 7T- 7 ] is usually negligible, it turns out that all terms in [r\ hm^\ are essential for the correct 
computation of the semiclassical equations of motion which are 

r = de/hdn — n x Q(n) 
M = eE + erx SB(r) - M.d5B/dr (30) 

where we have defined the vector B l = As shown in the following, the generaliza- 

tion Peierls substitution, not only is essential for the correct determination of the equations 
of motion, but also for the Bohr-Sommerfeld (BS) quantization condition. 



2. Bohr-Sommerfeld quantization 

To underline the relevance of the generalized Peierls substitution at the level of the 
semiclassical quantization of energy levels for an electron motion in an external uniform 



magnetic, we adapt the arguments of ref. [25| . For V = the equations of motion Eq. ( !30l 
become 

t = d(-^-\ and hn = eD ( ^-xb] (31) 
\ noil J \ noix J 

with D^ 1 = 1 + |B@. For convenience 5B = B is chosen in the ^-direction B =Bk, the 
energy reads e = Eo (tt)—M z (it) B. Consequently the orbits satisfies the conditions eq =const 
and 7r 2 =const. The semiclassical quantization of energy levels can be done according to the 
Bohr-Sommerfeld quantization rule 

K ± dR ± = 2n(n + 1/2) (32) 

where K± is the canonical pseudo-momentum in the plane perpendicular to the axis ir z = cte. 
The integration is taken over a period of the motion and n is a large integer. Now, it turns 
out to be convenient to choose the gauge A y = BX, A x = A z = 0. In this gauge, one 
has n z = P z = cte, and the usual covariant momentum Yl y = K y — j^-X. In this case 
the BS condition reads § H x dl\. y = ~\ eB (n + 1/2) . Assuming that the physically relevant 
variables are instead the covariant ones, writing thus BX = B(x — A x ) the generalized 
covariant momentum defined as 7r y = U y — ^-A x becomes 

n v = K y - ^-x (33) 
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which is formally the same relation as the one between the canonical variables, but now 
relating the new covariant generalized dynamical operators. This relation with the help of 
the equations of motion Eq. (13~T1) gives K y = n y + ^-x = thus K y is a constant of motion 
so that § K Y dY = K Y § dY = and Eq.([32]) becomes simply § K x dX = 2tt (n + 1 /2) . Now 
using the definition of the generalized momentum K x = tc x — ^-B and the differential of 
the canonical position dX = dx — dA x = ^| — dA x , the Bohr-Sommerfeld condition Eq. 
fl32|) becomes 

<j> n x dTT y = —2-neB (^n + - — -!- <j) A±dTT^j (34) 

where the integral is now taken along a closed trajectory T in the n space and j- § A±dnj_ = 
4>b is the Berry phase for the orbit T. Integration in Eq. (|34"|) defines the cross-sectional area 
§ n x dn y = S (e , K z ) of the orbit Y which is the intersection of the constant energy surface 
sq (it) =const and the plane tt z = K z =const. As shown in 25j Sq(e = Eq — M z B,tt z ) can 
then be written as 

„ . 2n\e\B f 1 , 1 rM z (n)dn\ , , 



here dn = ^Jdn x + diiy is an elementary length of the ir orbit. We have thus succeeded to 
deduce the required result Eq. (I33|) (first found by Roth [5| in a different way and without 
mentioning the Berry phase), as resulting from the generalized Peierls substitution in the 
BS condition. The importance in this expression of the Berry's phase for electrons in metal 



in connection to band degeneracy was later discussed by Mikitik and Sharlay Q. Eq. (35 
implicitly determines the energy levels E n (K z ). For instance for the case where the Fermi 
surface is an ellipsoid of revolution characterized by two effective masses, a transverse m± 
and a longitudinal mi one has 

to = tf(^ + ^-) (36) 



2m j_ 2mi 

and S(e,K z ) is a disc of radius square Therefore for M z =const we have the following 
generalized relation for the Landau levels: 

\e\Bh f 1 J mi \ h 2 K\ 
' n + - - (t> B - -r^M z 1 



mj_ \ 2 \e\ J 2mi 

which shows that both the magnetic moment and Berry's phase can influence the energy 
levels. A nice illustration of this result is provided by electrons in graphene with broken 
inversion symmetry jl71 ]. 
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B. Second order diagonalization 

As explained in section 2, the Hamiltonian diagonalization at this order requires the 
computation of the Berry connections and the covariant dynamical variables. This will thus 
be our first task. Note that, from now on, for the sake on simplicity, we will keep the 
notation introduced in the previous section, <5B = B. 



1. Computation of the Berry phases and dynamical variables to the second order 

Starting from the general expressions for the non projected Berry phases at the second 
order Eqs. I18H9B0I we have in the present case for A x = £A*da : 

A* = M*(V^K+^ 

^ 1-1 ( T, f 1 APlT7 _ , 1 



- i [,£ P • IV- < -A^ Kl e + -A?V Rl e + H.C. 



- l - [V-A% k V + A% 1 + H.C. } e klm eB m , X + hA, 



where now X = (R, K) and the same for Y. We aim now at writing the connection ^4 X 
in a more convenient form as an expansion in terms of the zero-order (actually first order 
in h) Af. To do so, we first start to express the crystal momentum Berry phase A as a 
function of the position Berry phase *4 R . The computation of A K involves a commutator 
with K + A^ that we compute first. Notice that the first order energy operator e i n Eq. 
(j28|) and the Berry connections at the zeroth order depend on II = K— |A as well as on 
Aq" 1 = eV R m A[AQ l , which allows us to write the commutator of K m + HAq" 1 with any such 
operator , e.g. [K m + hA^™, •] in the following way : 

[K m + Mf™, •] = eV Rm A t + hA^', 

Let us stress that this formula is not valid for the electric potential term appearing in the 
Berry phase formula since it depends exclusively on R. The action of \K m + HAq" 1 ,] on 
this term has thus to be computed independently. The additional contribution due to the 
electric field and the magnetic field to the momentum Berry phase is seen easily from the 
formula for the to be 

h 2 



Mo]" 1 • (A^V Rm V Rl V (R)) - | {V-Af x V + Ao} -V^B 
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Now the first term and the last term in the crystal momentum Berry phase can also be 
replaced as a function of the position Berry phase contribution. But at this point we need 
to go back to the technicality mentioned in the first section. While the final results do not 
depend on the way variables are symmetrized, one has to choose an initial symmetrization 
for the diagonalized energy at the zeroth order as well as for the Berry phase. For the 
diagonalized energy, we will choose to write it as a series of symmetric monomials in the 
momentum II. Concerning the Berry phase for the crystal momentum the natural choice 
is to symmetrize initially the variables such that at the lowest order in h one has A^ = 
\eV R m Ai (R) + \eA^V r^Ai (R). This will introduce some technical problem later, but 
it is the most simple choice for us. As a consequence, we obtain readily for the terms of 
interest : 



HA 



K 



h 2 



V + [A*] .V Y V + [A*™] + H.C.) 



+-eW Rm A l (V + [A*] -VyP + [At] + H.C.) 
so that can now gather all these terms and compute A Km as a function of A Rl : 



.h 2 



—i- 



, io}' 1 ■ (A^V Rm V Rl V (R)) - | {V-A^ x V + Af} V R „ 



B 



The potential Ai being developer to the first order in h. As expected, only the con- 
tribution due to the electric term is not rewritten as a function of A Rl . The quan- 



Mo] -1 • (A Rl V Rm V Rl V (R)) - f {V^Af x V + Af) .V Rri 
on the bands as being equal to : 



tity -if 



B 



can be computed 



h 2 



—%- 



Mop 1 . [A Rl V Rm V Rl V (R)J - - {V-JR x V + Af} .V Rm B 



.h 2 



= i- 



MN 



MN 



2 

if M 



N 



7 Rm V Rl V{R) h 2 

{P_A R x V + A*} MN .V Rm B HM^N 
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so that we have for the (M, N) component : 



AK m _ j 



+ 
1 

2 

+■ 
+ 



eV^A ( R+^aA A^+A^eV^A, U+^A* 



.h 2 {< l )V^ Rl V(R) h 2 



MN 



MN 



2 Sqm — £on 4 

[(eVi^^ (R) A Rl + A Rl eV Rm Ai (R))] 
^eVfl, V^A (R, P) (^?^f* + A Rfc ^) 



'MN 



MN 



.h 2 



•tf)v*-v*vw n 2 



MN 



^OM — e 0N 



-r{V-A*xV + A*} MN .V Rm B 



>MN 



where we defined the notation 5 MN = 1 — Smn, 

As explained before, the relevant variables for the physical expressions are rather the 
momentum variables II = K — A. Since at first order we already had A$ = eAf' x B we 
would expect a similar formula at the second order. However, since our method has been 
designed to deal with quantities symmetrized as functions of R and K, it appears that the 
momentum variable does not appear directly in our Hamiltonian. We thus postpone its 
introduction till the computation of the diagonalized Hamiltonian. 

We can now, as a second step, concentrate on the position Berry phase. We rewrite it as 



= hA* + ^ {A Rl W Rl A R + A^ l V K X + H.C) 



h 2 

f fl + JU^yMo]" 1 - 

- R + HAf, ^ (V-A Rn V + A Rr + H.C.) e nru .B v 



V- l l -A Rl V Rl e Q (X) + \a^V Ki s q (X) + H.C. 



We can also rewrite the third term (M, N) component as : 



i-i 



V- {^Vi^o (X) + l -AfV Kl e, (X) + H.C.y) ) 



1 \A Rl ) V Rl (e Q M + £qn) + [Aq 1 ) V Ki (som + £on) 

x V / MN V / MN ; 



if M ^ A 



if M = AT 



Using also, for quantities depending on II = K — A and R that \A Rl ) Vr + 

V /MP 

(.4^ ) = — eB. {(Aq) mp x Vn) at the lowest order in h we now give an expanded 
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formula for the position Berry phase at the second order : 

(A*) MN = (hA*) MN - g £ {eB. ((A*) Mp x V n ) (A*) pN + M ~ N} 

p 

h 2 (<') MN V Rl V (R) + eB. ((^) MiV x V n ) (s m + e m ) , 

PAr x Vn) (£op + £ojv) 

+y 1, (A )mp <W —3— 



& ^ { A » l ) MP VRy (R) + eB ^mp >< Vn) (bom + sop) 

d mp {A ) 

Fr\ ji T — K n o 



9 ^ 



2^ eoM-eoP — V~o, PJV 

"T Vn (WW x (Onn) - B *mn 
+4 (E «)mp («)p N x WU) -BW + H.C 

and its projection on Band N . 

V N (A%) = (A*) NN = (hA*) NN - | £ ( eB - («)np x Vn) (A*) pN } 

p 

4? (w- + KO KP w>™) 

4 (E (<*■ x Vn)) f^f^™ (A R ) PN - ff.C7. . 

+4 E «^)«p (W)p» >< (A R )™) " ((-4?) NP x (X)pp (A R ) PW ) B) 

P 

+H.C. 

For later purpose we will need also the modified Berry phase appearing in the dynamical 
Band variable r^. A direct application of our general set up of the previous section yields. 

A% = A% N + j (P N [AJ] .V Y V N [Af] + H.C.) 

R P R 

— A NN + — (A 0NN .X7r,A 0NN + A 0NN .VpA 0NN + H.C.) 
= A* N - (B. {(A*) NN x V n ) A* NN + H.C.) 

mately, the projected dynamical variables follow directly from our general method of 
;onalization. They are given by the canonical variables shifted by the modified projected 



Ultima 
diag 
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Berry phases : 



r N = R + A* 
k N 



K + A% = K + A% N 

n + a* = n + - 



dxB ^R+^,P+^<^ +H.C. 



Note that since — iV.R ro ^ [.,£o] 1 • [A^V^V (R) + H.C^j has no diagonal element 
'eVj^A (R) A R ' + <eV^A (R) 
+ yeV Rfc V^A, (R, P) P^ (^^J* + 



.4 



1 r 

2 



NN 



h 2 

8 
ft 2 



+ T eV« fc V* m A (R,P) (p w 



^H 



P 



A? 



4 



+ P 



iV 



./Hq 



P 



N 



s\ 



Band Hamiltonian 



We now turn to the Problem of deriving the diagonalized Hamiltonian. As explained in 
section II, the effective diagonal Hamiltonian for the Nth band can be written (see Eq. (pTj) 

) 



H N (X) 



e 0N (x) + ^Pn 



~P s v 



so (x),^ 
so (x),i Ki 



A Pl - 



so (x),i p < 



£o (*),A Pl 



A 



Ri 



,A R < 



V 4 



A Pl 
h 



e (x),A Rl ] A p > - [e (x) ,^]] A R \A Pl ] A R >] - - (e (x)> (37) 



Where P^r is the projection on this band. As explained in 21[ this Hamiltonian can be 
rewritten in a more enlightening way as : 

H N = e (tt) + M [e ] + [M [e ]] + V (r) - | (e (x)) 

Where the magnetization operator M acts on the following way : 

M [A] = -^V N I [a, - [a, A Kl ~\ A Rl } + E.G. 

As consequence, the Hamiltonian is given by a series of magnetization terms acting on the 
"classical" Hamiltonian, that is the Hamiltonian obtained for classical commuting dynamical 
variables. We will more give detailed formula for the various terms as functions of the first 
order Berry phases, but we turn first to the problem of choice of variables in writing the 
energy operator. 
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a. Momentum variables. Second order generalized Peierls substitution For the mo- 
ment, the pair of canonical variables implied by our method are the "non physical one" 
x = (r N ,k N ) 

r N = R + A% 



K + A% = K + A 



NN 



One aims at replacing k^ by a "physical" momentum variable which reduces to II = K— |A 
at the lowest order. To do so assume the electric potential is set to since it plays no role 
here, and rewrite the quantity Eon (x) that arises in the diagonalized Hamiltonian as : 

e 0N (x) = e 0N (k N - A (r N )) 

so that one would be tempted to choose — A(rjv) as the right momentum variable. 
However a computation of this quantity yields : 

k N - eA ( rjv ) = K + A* - eA (R + A%) 

= K - A (R) + l -e [vAt (R) A* + A R N l VA l (R) 

+y eV Rk VA l (R, P) V N (A* l A« k + A* k A% 1 ) 



+^eVV Rk A l (R,P) (V N 



V Rk A(R)Af + A%V Rk A(R) 



X 1 ] V N [A$ k ] +V n [A* k ] V N 



^eV^V Ri A(R,P) (V N 

h 



A* 1 V N A* k 



+ T 



N 



J[Rk 

✓-in 



N 



n + e<xB(R)+e- V N A% k (A* x V Rk B (R)) + B.C. 



+■ 



It 

8 



eV Rk VA l (R, P) V N yA l A k + A k A 1 J 



While the first three terms are physical since they involve the momentum II and the magnetic 
field, the last one is problematic since it will involve in the Hamiltonian some symmetric 
combinations of the vector potential of the kind V Rk Ai + V Rl Ak- The appearance of such 
terms is non physical but has nothing surprising here since we ordered all our expressions in 
R and K, not in R and II. As a consequence, we cannot expect in one individual expression 
to have only contributions of the momentum and magnetic field operators. An other way 
to understand this problem is to remember our choice of symmetrization for the crystal 
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momentum Berry phase that was easy to handle but does not fit when aiming at rewriting 
the Hamiltonian in terms of II. Actually having let the derivatives of the electromagnetic 
potential half on the left and half on the right automatically induces commutators terms 
proportional to Ai+V R,Ak. However, we know that the overall process of diagonalization 
of the initial Hamiltonian has to depend only on II, R, B (R) so that we know that the 
symmetric contributions V R k Ai + V Rl A k to the Hamiltonian have to cancel. 

This assertion will be checked explicitly for the above problematic 
term whose contribution (developing e 0A r (x) to the required order) is 
f eV Rk V Rn A t {R,P)V N (Ao l Ao k + A* k A*^) Vn n £ iv (x), when diagonalizing the Hamil- 
tonian. Actually we will find a counterpart to it. However, it is unnecessary in general 
to check the cancellation. These have to vanish by construction, so that we can discard 
them directly. Their local appearance in individual terms will be discussed briefly while 
computing the Bracket term for (x). 

As a consequence of this discussion, our relevant variable for the momentum in the 
Hamiltonian will be : 

7T = n + eA% x B (R) + \v N [A^ k ] {A* x V Rfc B (R)) + H.C.] 

and in the sequel, the relevant couple of canonical variables will be : 

(x) = (r N ,ir N ) = (r,?r) 

(the index N being understood). 

b. Magnetization terms At this order of approximation, and using our expressions for 
the Berry phases, one can directly develop the expression for M [e ] to obtain the following 
decomposition : 

M [e ] = -\(jjl (x) .B (x) + B (x) ./x (x)) - ^ (x) .V* fc B (x) 
+ Darwin term 



where the Darwin term refers to an expression which reveals to be of the same kind as the 
Darwin term in the Dirac Hamiltonian as we will see soon. 
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In this expression we have defined 



ih 2 



so (*),A Rl 
e (x) , A* 



A n " 1 1 



Jmn 



%h 2 



£o (X) ' ^ nm { X P ^ R )Jm, + - [ £ o (x) , { (V-A* x TVff),}] WM A% 



k 

OMN 



These expressions are directly obtained by expanding the magnetization term M [eq] = 



A Kl 



£o,A Kl A Rl | and computing the commutators that appear when 
letting the magnetic field half on the right and half of the left of the relevant expressions. 
Note at this point that, doing so, some contributions involving V R k A[ + V R t A k appear that 
cancel the corresponding contribution arising from the second order in H development of 
Eqn (x) as announced before. 

More precisely, the derivation of \x (x) and fi k (x) is as follows : starting with the term 
\V N { [e (x) , A R '] A p < - [e (x) , A p >] A R ' } we have to put the magnetic field half to the 
left and half to the right in the following way 



i 

T 



eV N { [e (x) , A Rl ] A Kl - [e (x) , A Kl ] A Rl } + H.C. 
= l -eV N { [e (x) ,A Rl ] {A^ViK + V t A n A R ^ - [e (x) , (i^V^„ + V t A n A R ^] A Rl ] + H.C. 

Letting the gradient of the potential to the left or to the right yields the required contribution 
for n (x) and fl k (x) is obtained by computing the commutators involving the gauge field. 
These terms are (we skip Vn for convenience and introduced a \ to take into account the 
integral over a) : 



16 



( [e (x) , A Rl j [v,A» - V n A h A Rn j ) - ^ [Vi A n , [e (x) , A R " 



A R ' 



16 

h 

eE 

16 



A*" [e (x) , V, A»] A* + [eo (x) , V, A,] A Rn A 



Ri 



+ H.C. 



nlk 



([e (x) Vn^V^flfc) - ^ [V,A», [^o (x) 



( [A R \ [e (x) , V,AJ + [e (x) , V^ n ] (i^i^ + i^i^)) + H.C. 

ViA n , [e (x),i R "j 



16 
ft 

16 



ee"' fe ([e (x) Vn^V^fl*) - ± ( 

' r A Rn , [e (x) , V/Aj] + [e (x) , V,A»] [i^, i Ri 
[e (x) , V,A»] + A Rl A Rn Yj + if.C. 



"I6 e 



+ 
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the projection of the last term 

[e (x) , V,A»] V N (A R "A Rl + A R 'A R ") + B.C. 
= -UviV m A n V N [A R -A Rl + A Rl A R ^j V Um e 0N (x) 

gives a contribution that cancel the one we have neglected in the momentum term as an- 
nounced before. Moreover 



-y^e[e (x) , VjAJ 



i*",^] +if.C. = - 



1 

—e [e„ (x) , ViA n ] , 
16 



A Rn , A Rl 



and this term does not contribute to the diagonal part at our order of approximation. 
The last term we are left to compute is then : 

~YQ e£nlk (h (x) Vn ^" v ^) - (iV [ V ' A - h (x) ^ Rn ] ARl \ + HC ) 



16 



i^Jeo (x),V,A,]l A R '+H.C 



16 



.raZfc 



s (x) , Vn p i R "l V^fifc) - [v,4,, [[e (x) ,A R -} ,A Rl 



e 

16 



A R \[e (x),V,Ajl,i*< 



16 

+ ^eV p V^„,Vn p 



= ~^ nlk ( [eo (x) , Vn^V^fl*) + A eVp V^„, V n , [ [e (x) , i*"] , A 

~e [A Rn ,[e (x),V^ n ],i^ 
= ~^ nlk ( [eo (x) , A*] Vn p A Rn V Rp B k ) 

,A Rl + e (x),A Rl 
A Rl ,[e (x),V ; A n 



R, 



Rn 



e 
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+ 



A R \[e (x),V,A, 
i*",^ ,[e (x),V,i4 



£ (x) , -4 



+ 



A Rl ,A R "^ ,e (x)]) 



+ 



,-4 



where we used twice the Jacobi in the last identity to rewrite 



^eVpVjA^Vn. 



e (X 



and — 

lb 



A R "\e (x),V,41,^ 



While iso- 



lating the contribution in the magnetic field, that is proportional to V R k Ai — W^A^ and 
introducing Vn only the first term in the last identity remains, which yields the first term 
in the magnetization jlf. The contribution V R k Ai + V_R,A fc can be disgarded as discussed 
in the next subsection. 
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The 



two 



last 



contributions 



for 



/V 



arise 



from 



the 



term 



— T {^-A? x 'P+ J ^} mn derived before in the expansion of A Kl . Inserted in 



- 2 eV N 



£o (x) 



,A R < 



A 



K, 



e Q (x) , A K[ A Rl \ yields directly the terms announced. 



Now, concerning the Darwin contribution, it is given by : 



Darwin term = — 
4 



. M 



h 2 



+-[e (x),A Rm ] 



MN 
R, 



V Rm V Rl V(R) 



NM 



Its 



derivation is straightforward and 



/ MN 

corresponds to 



isolate 



in 



\eV N | e (x) , A Rl A Kl - e (x) , A Kl A Rl } the contribution of A K > involving the 
potential term (see above). At our order of approximation, the ^4 X terms reduce to A^~. 

The last important magnetization contribution in the Hamiltonian is due to the squared 
action of the magnetization operator M. This is of second order in field and can be written 



M[M[e )} 



ih 



V 



N 



v 4 



v 4 



R, 



E , Aq 

V- 



$ A Rl 

En, J\. q 



An 



J^n 



Ri 



En, A 



A 1 



A* 



A 



Rn 



to 



, V-A Rk 



V-A R '),V-AZ m V-A^£ klp B p E mnq B, 



\Rr, 



A 



+ H.C. 



Rn 



c. Computation of —h{en(x)) The form for the diagonalized Hamiltonian given in 
(first equat) is taken from [21J. As explained in this paper the computation of the energy 
and in particular of the term — | (^o (x)) has to be performed by assuming a certain kind 
of symmetrization for the various expressions, that is a way to order the powers of R and 
K in expressions as En. Its value is in fact the track of the initial choice of symmetrization 
in the diagonalization process. Let us remark first, that since at the lowest order (that is 
zeroth order in h), only in (x) mixes R and K, the term — | (e (x)) will reduce to — | (in (x)). 
Second, and more importantly, let us stress that the final Hamiltonian is independent of this 
ordering, but this last one is necessary to give a precise meaning to each expression. The 
bracket term — | (in (x)) has in fact to counterbalance different choices of symmetrization 
and provide ultimately the same diagonalized Hamiltonian. 
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Concerning the choice of symmetrization between the canonical variables R and K a 
symmetric form in R and K is in general chosen, such as putting the powers of K half 
on the left and half to the right. An other choice is a complete symmetrization in R and 
K, that is permuting this to variables in all monomials, all permutations being equally 
weighted. However, both this type of symmetrization does not fit here since the relevant 
variable is 11= K — eA (R). We thus have rather to consider a symmetrization in the gauge 
invariant variable II. We will compute — | (eq (x)) in two cases that might be relevant for 
the applications, depending on the problem at stake. 

Assume first for Eq (x) a completely symmetrized form in the powers of components of II. 
That is, we will consider i (x) written as a series of terms of the kind v|> 111213 •••n ii n i2 II ia ... 
symmetrized in the indices ij. This form corresponds in fact to the most general one, not 
taking into account any spatial symmetry for the problem 

The (.) operation defined in [2l| has been recalled in the first section and can be computed 
explicitly for Eq, in the following way : 

Start with ^/ llt2l:i ^ ,ln Ili 1 Ili 2 Ili s ...Ili n , v|/* ll2 *3--- being a completely symmetrized tensor. 

Applying the rule given in the first section, the contribution of this term to — | (io (x)) 
is obtained by deriving (that is removing) by some momentum components IIj fc 11^ and 
inserting V ' A lk . One gets for this contribution : 

'4 Z V "" [ll,:-Hn : Vn r VA\ li . r ..U, i ...U„ -U^.M^.M,, : V„. . I'M I, ; . , ...1 1,., 

kl 

where the hat recalls that the variable is omitted from the series. Let move the gradient of 
the potential half on the left and half on the right to gather them by pair. We thus have : 



kl 



h 



-ie 



8 ^ 



l 2 ...lr, 



n ife+1 ...n ir ..n i 



kl 



h 



+ie-J2^ ili2 '" in 



ni 1 ...iii fc ...nj z _ 1 



^Ri,A ik - Vr-A 



n, +1 ...iT„ 



kl 



h2 V- w 



k<p<l 



h 2 v-7 



V Rl A lk +V Rl A ik 



n il ...n ife ...n ip ...n ii ...n i 



n il ...n ijk ...n ip ...n ij ...n 1 



k<p<l 



35 



the last equality is obtained since \& lli2 - ln is completely symmetric so that the contribution 
of Vb A %k — Vb A %k vanishes. 
As a consequence one has : 

— (e (x)) = —h (B k e jih ) Vn, Vn/o + contributions proportional to V CV j A i + ViA j ) 
2 8 

As a consequence, with our initial choice of symmetrization, — ~ {in (x)) contains only 
symmetric contributions. However, as we explained in the previous section, since by con- 
struction ultimately e (X) has to depend only on II and B (since B appears as the com- 
mutators of the components of II), this symmetric contribution has to find its counterpart 
in our previous development of the energy operator. This counterpart is automatically the 
contributions we discarded before. 

However, we can at this point explain a bit more how this works. Roughly, in the 
diagonalization process, the Bracket of the energy enters in two ways in our method since 
basically the trick to obtain the diagonalized Hamiltonian is to add and subtract (in (x)) 
(see [21], or more explicitly the differential equation in [14J]) to two different quantities that 
do not contain any symmetric term . 

The first way comes from developing 

£q (Xff + da + ^X Q ) — U a (X Q ) HQ (Xq) U^ (Xq) 

— U a (Xq +C 2q + g?Xq) Hq (K a+ da + dX Q ) [/+ (Xq+^q + g?Xq) — U a (X Q ) Hq (X a ) U a (X Q ) 

This term, yields all the relevant contributions to the diagonalized Hamiltonian we derived 
before. However, due to our initial choice of canonical variables R, K, the series expansion 
of {UH U + ) led us also to some contributions proportional to V (VjA 1 + VjA 7 ). They arise 
because the Bracket operation isolates non physical terms such VjA 1 . When we recombined 
them to get magnetic field contributions the V (VjA* + VjA 2 ) appeared. 
The second way comes from developing 

— [e a (X a+( fo, + dX Q ) — e a (Xq + ^q)] 

that will yield the contribution — (in (x)) including its symmetric contributions. But since 
U a (Xq) Hn (Xq) (Xq) and e a (Xq+^q) can be written in terms including only physical 
variables R and II (both the energy operator and diagonalization process are assumed to 
depend on these variables), the symmetric contributions of these two terms come only from 
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U a (X a+da + dX a ) H (X a+da + dX a ) U+ (X a+da + dX a ) - e a (X a+da + dX a ), and as such 
cancel by construction. 

As a check, one can see that, typically, the contribution arising from (UH U + ) were 
proportional to ViViVjU+Vi (V j A i + ViA j ) e (x) oc A^ViA^Vi (V j A i + WiA>) which is 
exactly the kind of symmetric terms we were left with while computing the magnetization 
term. But precisely these contributions are by construction equal to those produced by 
| (e (x)) (since UH U + = e (x)). As a consequence, — f {i (x)) acts in fact as a compen- 
sation term, and the overall terms in V (V ' jA 1 + ViA 3 ) is constrained to cancel whatever 
the form of the Hamiltonian or the Berry phases. 

As a consequence of this discussion, we can neglect the symmetric term and keep definitely 
for the choice of full symmetrization on the momentum components : 

-\ (e (x)> = 

However, the choice of symmetrizing all powers of the momentum II does not fit when 
the system presents some symmetries. The Basic example is the Dirac Hamiltonian 
satisfying rotational invariance, but one can also consider fermi surfaces having some 
ellipsoidal form. In that case a natural choice of symmetrization can be done with respect 
to some or several quadratic forms in the momentum of the kind U^A^Ilj (typically for 
the Dirac case A %0 = S lJ ). To inspect this case, we will thus compute — | (f (x)) when 
i (x) is written as a series expansion of products and powers of such quadratic forms 
HiA^Uj. To do so, we develop e (x) as a symmetrized series of monomial terms of the 
kind C (X) i U i (Q)U j (Y) j D. The (X) 1 and (Y) j being shortcut for 11, A*' and 
with A l i and some quadratic forms arising in the expansion. C, D, (Q) are arbitrary. 
Note at this point that implicitly, (X) 1 (Q) H,- (Y) 3 has to be considered as part of a 
sum of four symmetric terms obtained by permuting (X) 1 and H,- and (Y) 3 , that is 

\ [(xy Hi (Q) n, (Y) 3 + n, (xy (q) n 3 (y) 3 + (xy n, (q) (y) 3 h, + n, (xy (q) (y) 3 h/ . 

For the sake of simplicity we only keep (Xy Hi (Q) Uj (Y) 3 , the symmetrization being 
implicit and we skip the terms C and D that play no role in the sequel. (X) % Hi (Q) Uj (Y) 3 
contributes to — f (e (x)) by deriving (that is removing) H,- and inserting V 'jA 1 (the 
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sums over the indices i,j is understood): 

'-eh ((xy (Q) (vy - (xy (Q) (Yy) 

= \eh {{Xy V 3 A* (Q) (Y) j - (X) 1 (Q) (Yy) 

We can gather the potential terms by putting them half to the left of Q and half to the right 
of Q. We can thus write this term at order h 2 (which is enough for us here given our order 
of approximation) 

\eh ((xy (q) (y) j - (xy (Q) ( Y y) 
= jh (xy (Q) e^B k (Yy + |/i {xy ^ k B k (q) (Yy 

-leh 2 V t (VjA* + V^) (Xy V, (Q) (Yy + -eh 2 (X) 1 V, (Q) (Yy V, (V 3 A l + V^) 

8 4 
To the order h 2 , the terms f h(X) 1 (Q) e jik B k (Y) j + fh(Xye jik B k (Q) (Y)\ once summed 

over the symmetric monomials, involve all the possible permutations of e^ lk B k inside the 

series expansion of the second derivative of Eq. As a consequence, we can write : 

f% %P 

-- (e (x)) = — h (B k e jik ) V ni V n -£o + contributions proportional to V (V j A i + ViA j ) 
2 8 

As explained just before, it is understood that B k is inserted in a completely symmetric way 
inside the series expansion of V^Vie^o (that is cyclically permuted inside the series). 

Moreover, we can again discard the symmetric terms V (V ' jA l + ViA j ). 

We thus have with this choice of symmetrization : 

We will conclude this section by obtaining a developer form for l -^h (B k e jlk ) V^Vi^o- It 
is performed by going back to its series expansion. Let us assume in a first time that Eq is a 
function of II 2 only (rotational invariance), so that we can choose a natural symmetrization 
to express i (x) as a power series of n 2 . 

The sum of monomials of the kind 

IP ■ ■ IP 

-he (xy (Q) E^ k B k (Yy d + -he (xy ^ k B k (q) (Yy d 

can be written (we now reintroduce our implicit permutations that let B k to be half on the 
left and half on the right of (Xf and (Y) j ) : 

IP IP 

-hCU t E^ k B k (Q) UjD + -ftCTL, (Q) E^ k B k U,D 

IP IP 

+—e iik B k fiCU i (Q) TljD + —hCTl t (Q) Tl 3 E^ k B k D 



38 



Letting move the magnetic field half to the left and half to the right yields to the order h 2 
yields : 

^hCU t e^ k B k (Q) UjD + ^hCU, (Q) e jik B k UjD 
+^ ifc £ fc /iCTI i (Q) UjD + ^CU t (Q) UjS jik B k D 

IP IP 

= e ltk B k -h [CUi (Q) UjD + CU { (Q) UjD] + -h [CIU (Q) U 3 D + CTT (Q) UjD] e jik B k 
8 8 

+ (V x B).n-h 2 C{Q)D 

8 

Concentrate on the first term. One has 

IP IP 

e ]lk B k -h [C^ (Q) U,D + CU { (Q) UjD} + -K [CUi (Q) U 3 D + CU t (Q) UjD] e jik B k 

o o 

7 P 

= e jik B k -h [CIT (Q) UjD - CU 3 (Q) U t D] 

o 

to the order h 2 . Now, 

7 P 

e jik B k -H [CTT (Q) UjD - CU 3 (Q) U t D] 

O 

7 P 

= jhC(Q)[U t ,Uj]De^ k B k 

IP IP 

+ jHCU t [(Q) , Uj] De^ k B k - -HCUj [(Q) , IT] De> ik B h 

Summing over the symmetric monomials of the series expansion of Eq , each of this commu- 
tator will yield a corresponding contribution to — | (i (x)). The first commutator will yield 
the contribution : 

^hC(Q) [Ui,Uj]De jik B k = ^h 2 C (Q) D (e ijk B k ) 2 

o o 

- ~h 2 {e^B k ) 2 A{V u ,) 2 e Q 
= e l h 2 (B 2 ) (V U2 ) 2 e 

the global | amounts for the fact that given our conventions the derivative with respect to 
IT are always on the left of the derivatives with respect to IT,-. This implies a \ factor each 
Vn^V^ term. We have also used the fact that given the rotational invariance (VnJ 2 and 
(VnJ 2 can be replaced by 4(Vn 2 ) 2 - Actually the commutator [IT, IT,] amounts for taking 
twice the derivative with respect to IT and IT,- (times e ljk B k ) and twice the derivative with 
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respect to Ilj and IT, yield each a contribution proportional to Vn 2 • On the other hand the 
second commutator will lead to the contribution : 

ie , „„ r „ „ ,„•!,„ 1 e 2 



+jhOn j [U i} (Q)\ De^ k B k -> -^ft 2 (v^Vpx.2 (V n *) e ) e^e*"^ 

-h 2 ((5 2 ) 2V ri .V ri . (VnO^o - B k V nk B m V Um 2 (V n 



e 2 
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with the implicit convention, recalled by the ', and implied by construction and the defini- 
tions of (Q), Ilj, Hj), that a factor II 2 that has been derived by Vn m will not be derived 
again by Vn fc (the same convention applying of course for V n; V' n ^). Similarly V^V^. 
means that a power of II 2 is not derived twice by Ilj. 

The reason for the global ~ factor is the same as before, the derivatives with respect to 
Ilj being on the left of the derivatives with respect to Ilj. A | factor has also arisen from 
the fact that the derivative with respect to Hi has to be taken between the derivatives with 
respect to Ilj and Ilj. Due to the symmetrization of the variables, it amounts for only one 
third of the derivative with respect to II/. 

Similarly, the last term will give : 

2 

~hCXl j [(Q),Th\De? ik B k - "-fi 2 (V' n , V Ui (V n ») io) e ijk B k £^ m B m 

,2 



^h 2 ((5 2 ) 2V rij V / n . (V n2 ) 2 e - B k V n B m V nm 2 (Vn?) e„ 



with the same convention as before for the derivatives. 
Ultimately, we compute similarly the contribution : 

{V xB).U e -h 2 C{Q)D -> i^ 2 (VxB).n4(V n2 )^o 

o z o 

= ^ 2 (VxB).n(V n2 ) 2 £ 
Gathering all the relevant terms yields in the end : 

fc 2 2 

- 2 <*o (x)> = e -h (B 2 ) (V n2 ) 2 4 + ^h 2 ((5 2 ) V n .V n . (V n2 ) 2 - B k V n B m V Um (Vn?) e, 
+|fr 2 (VxB).n(V n2 ) 2 4 



A direct computation shows that specia 



electromagnetic field, will be given in 221 ] . 



izing to the case of the Dirac Hamiltonian in an 



The more general case can now be treated in a very similar way. Coming 
back to the general form for the monomials, leads to consider again the monomials 
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-fh(Xye^ k B k (Q) (Y) j - fh(X) j (Q) e^ ik B k (Y)\ Here to alleviate the notation we as- 
sume again implicitly that B k is put half on the left and half on the right of (X) 1 and 
(Y) 3 ). the contribution to move them on the left of (X)* and the right of (Y) 3 leads to a 
similar contribution to — | (f (x)) as before : |ft 2 (V x B) .Vn (Vn 2 ) 2 £o Concerning the 
other terms, we recombine the monomial — ^h(X) % (Q) (Y) 3 £ ljk B k with the corresponding 
contribution +f h(Y) j (Q) (X) 1 e ljk B k (the e jlk B k being implicitly half on the left and half 
on the right of the expressions): 

IP ■ ■ IP 

-jh (XY (Q) (Y) 3 e i3k B k + -h (Y) 3 (Q) (Xf e i3k B k 



■ ijk B k 



(xy (Q) (y) 3 - (y) 3 (Q) (xy 

= ((Xy (Y) 3 - (Y) 3 (Xy) (Q) e i3k B k 

~(Y) j ,(Q)\ (xye i3k B k 



8 



(X)\(Q)\ (Yye^ k B k + -h 



As before, summing over the symmetric monomials of the series expansion of io , each of this 
commutator will yield a corresponding contribution in derivative of Sq. The first commutator 



m 



(Xy (Y) 3 — (X) 3 (Y) % is computed in the following way The rotational invariance does 
not exist anymore now, but (X)* and (Y) 3 are still functions of II. As such they will yield 
the contribution : 



le 
8~ 

2 8 



h(Q) [(xy (Yy - (xy (Yy] e i3k B k 

h 2 {e i3k B k ) (e lmn B n ) (V n J (V n J (V n J (V n ,) e 



Once again, the global | amounts for the fact that given our conventions the derivative with 
respect to IT are always on the left of the derivatives with respect to IT,-. 
Similarly, the second and third commutators will lead to the contribution : 



ip r ■ ~i ip r 

■-h(Xy [(Q) , (Y) 3 ] e* 3k B k + -h(Y) 3 [(Q) , (X)' 



Inm r> 
Dr, 



■ ijk B k 



- T / (V hi V^ [(VnJ . (Vn,)] e ) e^B k e 

~t/ ( Vn ^ [(VnJ • (VnJ] io ) £tJkBk£lnmB ™ 

with again the implicit convention, recalled by the ' and implied by construction and the 
definitions of (Q), (X)\ (Y) J ), that the derivatives V n and V n . are not applied on the same 
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quadratic term HiA^Hj in the series expansion of i . [(Vn„) • (Vn.,)] is a notation to recall 
that (Vn n ) (Vn-) ac t both on the same quadratic term, 

Gathering all the relevant terms yields ultimately for the general case : 



h 



a (io (x)> = e -n 2 {e ijk B k ) {e lmn B n ) (V n J (VnJ (V n J (V n J e 

+YA h2 ( V * V * [( Vn ") • ( Vn ^)] io ) ^ kB ke lnm B m + 6 -h 2 (V x B) .V n (V n ,) 2 e 

d. Final form for the Band Hamiltonian Gathering all the previous terms, and writing 
io (ir) as a completely symmetrized series in the powers of the momentum (in the lack of 
any a priory symmetry), we have ultimately the diagonalized Hamiltonian : 

H d = e (tt) + V (r) - ~ (fi (x) .B (x) + B (x) .ll (x)) - fi k (x) V Rk B (x) + B./LB 



h 2 



h 2 



An 



MN 



M 



+ -[e (x),<-] 



MN 



V Rm V Rl V(R) 
W Rm V Rl V(R) 



NAI 



(38) 



The double scalar product B./2.B meaning that the two index tensor Jx being contracted 
twice with B. 

e. Example : Darwin term for the Bloch electron We now consider more specifically 
what we have called the Darwin term. Its interpretation turns out to be more transparent 
when B = 0, therefore we keep here only the electrostatic potential. In this case, we have 



A 



R 



4 



li 2 



and Aq 1 = as shown in 



h[.,e )-\(v-^AfV R y(R) + H.C.yj ,R + A t 



161 ] but at the second order we have the following contributions 



and 



(-4") 



MN 



2 - rt « 



MN —- S 



£qm — £qn 



MN 



showing that the non diagonal part of ^4 K is non null. It will lead to a magnetization term 
: the Darwin term. We can write : 



M R U = (^) MN + ^{(AS-) Mp V Rl (A^ PN + M»N} + t l l V K [A °' 



MN 



Vp^(R) 



h 2 

+ Y ( A o*)mp 5pm 



A Rl 



PAT 



VaV(R) 



h 2 ' & 



MP 



Vp ; v(r; 



£op — £qn 



£om — ^OP 



£qm — £qn 
<W (A*) pn + H.C 
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and the projected variables : 

(A K ) N = (A R ) NN = (M*)„„ + | (,(*) Vh« K) 



PAT 



£ 2 

+ y 



«)N P {^) pN +{A^) Np (A^ pN )5 PN 



£op — £on 



Therefore the Hamiltonian is : 

H d = i (k) + V(r)+M[e } 

= i + V(r)-^V + {[[io,A^} [X,K m \] - [[e Q ,[X,K m }]A«™}}+H.C. 



with : 

which reduces to 
H N 



X=[.,e 



-1 



.(v^^V Rl V(R)+H.cX) 



h 2 

i + V (r) + - 



M 



NM 



•4? i x />' 

1 MN 



V^O > MN 

V R ^V Rl V(R) 



v Rr y Rl v(R) 



NM 



+ H.C. 



(39) 



£0M — ^OiV 

This is an interaction in second derivative of the potential, Darwin type. 

The effective in-bands Hamiltonian Eq. fl38l) is the desired result concerning one Bloch 
electron in an external electromagnetic potential and will be very useful for the latter com- 
putation of the effective Hamiltonian of several Bloch electrons in self-interaction. But at 
this point it necessary to compare if it is possible with previous approaches in particular 
with Blount's one. 

/. Comparison with previous results In a powerful series of paper, Blount 



4j(27| de- 
signed a method to derive the diagonalized Hamiltonian both for a Bloch electron in a 
constant magnetic field and the Dirac electron in an arbitrary magnetic field as series ex- 
pansion in powers of the fields. In both cases, at the first order in h our results coincide 
with his ones. At the second order, the direct comparison for the Bloch electron is difficult 
to do, since Blount's choice of variables (the canonical ones) differs from ours. However, 
the comparison in the particular case of the Dirac electron (2 bands Hamiltonian) can be 



performed and will be explained in detail in 



22j | . The result is that, despite some important 
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similarities, our results differ slightly From Blount's ones. Actually some mistake arise in 
his results due to his choice of variables, which induces wrong expressions for the Berry 
curvature, and as a consequence, in the Hamiltonian. As a consistency check we have shown 
in 22], that at order -k-, we recover the usual Foldy-Wouthuysen formula for diagonalized 
Hamiltonian of the Dirac electron in an electromagnetic field. This is apparently not the 
case in 27| where the coefficient of the squared magnetic field is apparently incorrect. 



IV. INTERACTING BLOCH ELECTRONS 



Having found the diagonalized Hamiltonian of a one particle system, we can now focus on 
the multiparticles Bloch electrons in interaction. In this context we must consider the inter- 
action mediated by the Coulomb interaction. Although this last one dominates the magnetic 
interactions, it is known that, for material whose electrons have a non vanishing magnetic 
moment, new effects such as interaction between the moments mediated by the magnetic 
field can have important effects. For this reason, we are now interested in investigating the 
diagonalization of the Hamiltonian of several electrons in self-interaction through the full 
electromagnetic field. However, for practical reasons and to deal with tractable formula, we 
will restrict ourselves to the case of relatively small interaction. 



A. Derivation of the microscopic Hamiltonian 

To start , we will derive the microscopic classical Hamiltonian for this system, before 
going to the quantum version. Note that the electromagnetic field will always be considered 
as classical and only the particles will be treated at the quantum level. Let us introduce the 
non relativistic Lagrangian of P particles self-interacting through the electromagnetic field 
and moving in a periodic potential V p which is the only external potential: 

L = ^J2 mia) ( rW ) 2 - V p ( R(Q) ) - \ F ^ F>XV + [ d 3 x3.A- I d 3 xeA J 

a - J J 

Here is a classical electromagnetic field whose potential is (A,Ao), e are the particles 
charges. The current J is given by J (x) = J2 a e6 ( x - R(Q) ) R(a) and J o is the density of 
charges, J = J2 a e5 ( x _ R(a) ) ■ 



44 



Choosing a gauge for the electromagnetic field, for example the Coulomb gauge (see [28] 
), some computations lead to write the particles and field Hamiltonian as 
p 



H = ^P (Q) R (Q) - i^m (Q) (R (Q) ) 2 - J d 3 x3.A+ J 

a=l a 

-l(U + VA ) 2 + V p (R^) 



(f\rJ.A+ I t f 6 xeA J + n^ m + i (V x A) 2 



with il em the solenoidal part of the electromagnetic momentum defined by Il em = A satis- 
fying the constraint V.Il em = 28j- is the usual canonical particle momentum. Using 
also our definition for J, we are thus led to : 



' (pW-eAW(RW)) 2 i 



a=l 



2m( Q ) 



efxe 



A (x) Jo (x) + l -Iil m + I (V x A) 2 - (VA ) 2 



+V P (R^) 

^(pW-eAW (R (a) )) 2 + 1 



a=l 

+ 2 ' ! '' 



2m( Q ) 



d 3 xeA (x) ^ e5 ( x - R 



n 2 m + o d 3 x(V xA) 2 - (VA 



+ V p (R< tt ») 



(40) 



Recall that V^, (R^) stands for the periodic potential. Note also that, given that the particle 
density is a sum of delta function centered around the particles positions, the potential 
\ J d 3 xeA (x) £ Q ed (x - R (a) ) reduces to : 

\ J (x) £ rf (x - R<">) = £ 



However, we will not use this "solved form" right now and rather keep the electromagnetic 
field A (x) (see below). 

The previous expression is desired Hamiltonian for P classical electrons in interaction 
and in a periodic potential. It will be the starting point for the quantification and for the 
computation of the effective in-bands Hamiltonian to the second order in h. This order is 
required to determine the magnetization-magnetization interaction in the same way that the 
spin-spin interaction in the Breit Hamiltonian for Dirac particles 23|]. The consideration of 
interacting electrons, each of them being constrained to be in a single band, will require the 
Hamiltonian diagonalization and the projection on each respective band. To do so we just 
have to apply the method developer for the one particle scheme. 
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Remark, at this point that the Hamiltonian could be further simplified by using the 
Coulomb gauge condition which allows to replace directly the potential A . Actually, as a 
consequence of the Coulomb gauge constraint one has — V 2 Aq = J . However, for the sake 
of symmetry we aim at treating the electrostatic potential as an external field, as we do for 
the vector potential, and to replace it after transformation. We will explain later why this 
is in fact simpler and innocuous for the final result. 

B. Diagonalization process 

We will now diagonalize our Hamiltonian to the second order in h as before. The space 
of states acting is the tensor product of P copies of the individual spaces V = <8> a =i,...,pV Q , 
with the Hamiltonian Eq. fT40l) in which V p is the periodic potential, is the potential 
created by the other charges on charge a and v (R (q) ) = / d 3 xeA (x) e5 (x - R (a) ) the 
electrostatic potential involving the particle a. 

The Hamiltonian H is diagonalized at the second order in h in a straightforward way. 
Using our previous experience with the one electron, consider : 

Where are the individual particle diagonalization matrix found in the previous section 
to the first order in H. Note that this matrix do no act on the electromagnetic field, whose 
integration variable x is independent from the quantum mechanical operators Ft*"-*, p( a \ 
This is the usual procedure for this kind of system. We first diagonalize the part of the 
Hamiltonian which describes the particles dynamics and leave the electromagnetic field 
untouched (see for instance Feynman's book for this procedure |29|). This is essential in 
order to have each particle living in a particular energy band. It means that we assume here 
an adiabatic process in which the electromagnetic interaction does not cause interband jump. 
If we were eliminating straight the electromagnetic field in the initial Hamiltonian Eq. ( 1401) 
with the help of the Maxwell equations, we would of course get an equivalent Hamiltonian 
for particles only, but it would automatically mix all the energy bands. Then the possibility 
to assume an adiabatic process and to project each individual particles Hamiltonian on a 
specific band would be lost. 

Therefore we will choose the first route and after the diagonalization one will solve the 
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Maxwell equations which will determine the exact form of the particle currents. These last 
ones will then be further injected into the Hamiltonian to get the final in-bands effective 
Hamiltonian of the particles. Let us show how this procedure works. 

From Eqs. ( 1551) ( 1401) the Hamiltonian in the diagonal representation of the particles part 
is thus directly : 

h n = J2 4% + E v ( r(Q) ) - \ E - B(a) + B(a) * w) 

-jl k (x) .V Kfc B (a) (x (a) ) + B (a) ./2.B (a) 

-Af(«) 



+ - 



M( a )N( a ) 



+tE 



4 "* ( x ) ' -A) 



(a) 



A R i 
-H 



(a) 



(a) 



a) 



(rW) 



(a) 



\2 



+1 J d 3 x Tl 2 em + \J d 3 x (V x A) 2 - (VA)) 

The notations of the previous section have been kept, adding only the superscript (a) to 
label the particles in Eq. (1551) . For example, x^"^ = (r( a ),k( Q )) is the couple of dynamical 
variables for the particle (a), x*") = (tl^ + A^], + A^IJ^ and the Berry phases 
^^(a) , ^Tvf") involved i n these operators have been computed in the single particle section. 
Notice also that N is a shorthand for the multi-index iV = (iV^) recalling that the 

Hamiltonian describes P particles respectively in the bands N^\..N^ P \ Ultimately is 
the field acting on particle a. 



C. Effective Hamiltonian for P particles and particles currents: 

Our aim is now to replace the electromagnetic field in the Hamiltonian as a function of 
the dynamical variables of the particles system. We first do so by solving the Hamiltonian 
equations for the electromagnetic field in the limit of relatively weak interactions. 



1. Dynamics for the Electromagnetic Field 



To get rid of the electromagnetic field and writing an interaction Hamiltonian for the 
particles, we first focus on the vector potential and start with the Hamiltonian equation for 
the electromagnetic field in the Coulomb gauge |28| : 
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A(x,t) = n em 

KmM = ~Jd 3 y 



5 ij 5 3 (x - y) + 



d 2 



1 



dx l dxi 4ir |x — y| 



- I d 3 y 



H m at, ^Im (X, t) 

5 ll 5 3 (x - y) + 
5 lJ 5 3 (x - y) + 



d 2 



d 6 y 



dx i dx^ 4ir |x — y| 

d 2 1 



V x VA j (y, t) 

V x VA j (y, t) 

5 



(x, t) 



H, 



mat 



dx l dxi 47r |x — y| 

where H mat denotes the matter part (that is excluding the free part of the electromagnetic 
field Hamiltonian). Using the gauge constraints : V.II em = V.A = 0. we are led to : 



□A 



((x,t)) = - J 



d 3 y 



S ij S 3 (x - y) + 



d 2 



1 



dx l dx'j 4tt |x — y| 



8Ai(y,t 



-H, 



mat 



where t - ) stands for the functional derivative with respect to A % (y, t) 



2. Explicit expression for the current 

We aim now at writing a more explicit formula for the effective current 
H mat . To do so, we decompose H mat as : 



mat 



5A?(x,t) J 

EC (* < ° > )+E»(r ( ° ) ) 

a a 

-\Y.{v ( x(a) ) - B(a) ( x(a) ) + B(a) ( x(a) ) ^ w) - E ^ w • v ^ B(a) ( x(a) ) 

+B (Q) ./2.B (a) + 



where stands for the Darwin term : 

E Uf) (^) 

'[4 a) (x),<™ 



£>(«) = ^_ 

4 



i-R 



(a) 



I V R ( a )V R (c) 
iV(«)M(«) V / M(a)iv («) m ^ 



(rW) 



The computation of j J (y, t) = — 5A j S ^ t ) H mat has to be performed carefully since the depen- 
dence in the gauge field in H mat is intricate. We aim at computing it as a function of the 
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"physical" variables r^ a \ and not as a function of the initial microscopic canonical vari- 
ables. However, the physical variables r( a \ tc^ depend on the value of the fields evaluated 
at the microscopic positions and momentum. 



r («) = R w + ^ (a) ^R(«),nw) 

n (a) = n (a) +eA f^ ( R W R^) X B (RW) 



+e-V N 



R 



A k 



(a) 



Af a) (RW.nM) x V R (.)B (RW)+EC.) 



to turn this difficulty we reintroduce the variables r*"), 7r( a ) recursively in our Berry phases. 
Starting with Y^ a \ at our order of approximation and using R/ Q ) = r^— A^ a) (r^, n^) 



= rM+1 RM ( r W,7r< a >) 

(we have neglected the terms in V r («) that do not appear at order ft 2 and skipped the 
index N for the sake of simplicity). We do the same thing for the momentum by start- 
ing to write it as a function of k.( a \ v( a \ Given the definition of the momentum, we 
explained before that it was not given by k( a ) — eA (r^), but rather k( a ) — eA (r^) — 
^eV Rk VAi (R, P) V N (A* l Ao k + -AJMJ 1 ) . It was explained that this variable had a phys- 
ical meaning (the subtracted term was there to compensate a non gauge invariant term). 
This last term cancels one term in the Berry phase for K^. As a consequence, 

n (a) = k («) _ eA ( r (a) j = K (a) + ^ R (a) > n (a)^ _ gA ( r («)) 

= K (a) - eA (r (a) ) + ^e [vA (R (a) ) A# (RW,n (a) ) + A% (R^ a \n^ VA t (RW) 
ft 2 



+-eVV^ (R, P) Pjv [< fc ] + P W [< fc ] V N [a* 1 ]) 



= K( Q ) - eA (r^) + -e VA, (R (a) ) 2^ (r^, vr^) + T l (r^, tt^) VA, (R (q) ) 

n 2 



— eVVj^, (R, P) (V N [A* 1 ] V n [A^ k ] + V N [A^ k 



V 



N 



A Rl 



K H - e A (r (a) ) + -e 
+ft 2 eVV^ (r^) A lk 



VA t (R^) A Rl (r^,^) + A* (r<- a \n™) VA t (R (a) ) 
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Actually, as seen on their definition, the variables ir^ and depend on A (r^). Rather 
than solving these circular equations where the variables of interest appear on both side as 
functions of the canonical variables R/ a \ we rather unsolved the dependence of the 
variables in A and keep this circularity to compute the derivatives of interest, since it will 
allow to express all the results as functions of the variables r^ a \ 7r^ a \ 

Keeping this in mind, we start by considering a simple example that will allow to find 
some general rules for the computation of the functional derivatives with respect to the 
gauge potential. Let G be a function depending on a variable x (in our context r( a \ ir^) 
depending itself on a function a of a; (here the gauge field) through the relation x — f (a (x)). 

We aim at computing first j^^G (x) where y is an independent variable (i.e. the space 
parameter in our context). 

Start first with = {a (x)) : 

6 -*='/(«(*)) = ^ («(*)) 



Sa(y) Sa(y) 5a (y)' 

It would be wrong to conclude that reduces to 5 (x — y) because y is an independent 
variables but x is not and depends on a. One rather has to consider that a (x) is an infinite 
series of composition a(x) = a (/ (a (x))) — f (a (/ (a (x)))) — f (a(f (a (/ (a (#)))))) and so 
on. A slight functional variation of a propagates along all the series and has to be taken 
into account to write symbolically the infinite series : 

Sa (x) = S[a(f(a(f(a(f(a(f(a...] =8a(f(a(f(a(f(a(f(a... + a(f(6a(f(a(f(a(f(a... 

+a(fW(Sa(f(a(f(a... + a(f(a(f(a(f(5a(f(a... + a(f(a(f(a(f(a(f(5a... + ... 

where the inserted 5 acts solely on the a directly on its right. The variation Sa (x) can be 
rewritten as : 

5a (x) = 5a (x) + a(f(5a (x)) + a(f(a(f(5a (x)))) + a(f(a(f(a(f(5a (x))))) + ... 

where now x can be seen as a frozen variable (that is not depending on a) everywhere in 
the right hand side. A direct application of the chain rule yields directly : 

= 6(x-y)[l + a' (x) f (a (x)) + (a' (x) f (a (x))) 2 + .. 

so that : 



5 8a(x 



5a (y) 5a (y) 



x = f(a (x)) Pf( = 5(x-y)f (a (x)) 1 + a' (x) f (a (x)) + (a' (x) f (a (x))f + 
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and ultimately : 

S 



5a (y) 



G (x) = S (x - y) a (x) f (a (x)) 1 + a' (x) f (a (x)) + (a' (x) f (a (x))) 2 + 



Of course, this expansion is formal, and assumed to converge. Moreover, in the sequel the 
sum will always be truncated to a finite order. 

We can now translate this results in our context to compute the derivative — t ) H ma f 
since it involves only a generalization to several variables. The formula are a bit more in- 
volved since the circularity depends on three variables. Actually depends on through 
the Berry phase, depends on A, that depends on r^ a \ products like a' (x) f (a (x)) will 
be now replaced by products of three types of derivatives, — VA. How- 

ever, a simplification arises here. Actually, since the beginning, we have assumed that the 
magnetic part of the interaction between the particles is relatively weak with respect to the 
electrostatic potential. We will thus assume that the current is relatively weak, and the 
expansion will be performed only at the first order in the field. As a consequence, at this 
order of approximation, it is useful for the sequel to note that the composition series for 



<5A*(x,t) 



reduces to 



8 k + V t A k (rW) 



dr 



dA*(r(«)) 



Introducing the needed indices, taking into account this approximation and going to 
the second order in K we are led to the following results for the various terms involved in 



5 TT 

"<L4J(x,t)- nma * • 



-E 

a 

= -E 

a 

and 

-E 



6 ' 45 («) 



£)(«) 



6A i (x, t) 



d 



dn 



(a) 



X,t 



4% 



dir 



(a) 



dA k (rW 



dA i (r(«)) 



5(, 



5 A 1 (x,t) 



(a) 



dv (r^) dr 
Q r ^ dA k (r( Q )) 



5 k + V t A k (r^) 



dr\ 



(a) 



dA i (r( Q )) 



5 (x - r (a) ) 
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and also 



5A i (x, t 
+B (o) ./LB (a) 



-\ E ( x(a) ) - B(a) ( x(a) ) + B(a) ( x(Q) ) w) - w - v ^ B(a) ( x(a) ) 



IE 



(a) 



d 



+ 



a) 



d 



dA* (rH) 97r (") ( X) t ) dA* (rW) 0> 



(/x (x^)+/x fe (xJ.Vj .BW (x<°>) 



x5 (x-r (a) ) +H.C. 



1 V- 



(/. (xW)xV) Uf+V^ (rW) 



*7 



(a) 



5 (x - r (a) ) 



dA* (r(°)) 

+ [/t fc (x^) xV Rfc V<5 (x - r (a) )] .] + if.C. 
! [b (x (a) ) ./I (x (a) ) xV(5(x- R (a) ) - V5 (x - R (a) ) x /I (x (a) ) .B (x (a) ) 



These expressions involve some derivatives that are computed as follows, as directly implied 



by our previous remarks. d ^ r ( a )^ > rf^r(°)) cover ^ ne dependence of the transformed dy- 
namical variables in A 1 (r^)). Recall that they have to be computed such that in A 1 (r^), 
the is frozen, so that A 1 (r^) is seen as an ordinary, or independent variable. 

Starting with dAt ^ r(a) ^ the dependence in the field comes from the Berry phase of r\ , 
A R< i ) which is a function of the two variables 7r^ and r^ a \ We can thus write : 



» 



dr 



(a) 



OA 



+ 



OA 



R 



(a) 



drt ] 



+ 



OA 



R 



(a) 



dA i (rW) \ dA i (r( Q )) cL4* (r( a )) dA* (r( Q )) 



,(«) 



the second term in the right hand side is of order h 3 since 9A , is of order h 2 and — dr , n l 

to dr ( ™> dA l (r(«)) 

is of order h. At order h 2 one thus has : 



dr 



(a) 



;R 



(a) 



OA 1 dir., 



(a) 



&4 



(a) 



+ 



^(rW) faM dA^v^)) dA i (r^) 



,(«) 



Practically, the partial derivative 9 ^ r ' (a) ^ w hh respect to the field is obtained by decom- 



posing : 



OA 1 



~R (a) 

dA 1 d 



- (AT'.V^Af (rW.TrW) 



&4* (r( Q )) 



R (a) ( ) 

rM-(V r (,xV lW )4° (r<°>, 
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(a) 



and d^f/ r '( Q) \ is carried by the following term computed previously in the one particle case : 



^r = 4E{eB.(K w ) WF xv nw )K"0 PK } 

-4 (B. ((A RW ) m x V n-..) -4*li + H.C) 



F 

+H.C. 



8 , 



dAfo) i s a two tensor whose components (i, I) is given by : ^j^pj — e ljk X i [a 'Vk 



,(«) 

, i a; 

NP \ J PP 



A* (a) ) x(A* (a) ) (jtf ) ) \5 PN + H.C. 

PN/ j 



and the gradient acts on 5 (x — r^). 

Having given some expanded formula for - . j . , ? it thus remains to compute - n , n . , N . 
To do so, recall again that the variable is "frozen" with respect to A, which is thus 
thought as an independent variable. We use also our previous result on the Berry phase : 

^ = - eA (rW) + \e [v m ^ (RW) A* + ~A Rl V m A t (RW)' 
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Define also A mi = e-ferV r ( a) V( a )A r (r^) + e k n>i V p (a) Vr«) so that one has : 



d~A 



ecL4*(r( a )) 1 r ™ r - v y cL4 l (r(°)) 

» 



fA 



m.i 



= ^ _ A^Vra, - V< a) A n ( r W) + A 

* r^' V / (r( a )) 



faW ^'(rW) a r (°) ^4 Q) <*A* (r(«)) (r(«)) 



^ (r(«)) + a^wy 1 ~ ' 4 ' 



to our order of approximation. 

We can rewrite the last equation as : 



d-Ku 



(a) 



dA i (r( a )) 



5™-eV r ( Q) A n (r< a >) 



dA " 



97T. 



(a) 



e - 1 R ^V (a) - V {a) A n ( r W) 



(r(«)) 



+ A 



which is solved to the second order in h and first order in field by : 



1 dirtf 
edA i (r( a )) 



I'm 1* f7i * ' 



&4 



(a) 



m * m 



(r(«)) 



9A ' 



= 8? - 1<V (a , - V (a) A n ( r W) ' 



ft 



(a) 



(r(«)) 



(a) 



+eV r(Q) A" ( r W) 



9 y4 



(a) 



dir. 



(a) 



"A 1 V fa) + A mA 



(a) 



the gradient appearing at the end of the right hand side has to be understood as acting on 
the delta function 5 (x — r^) in the expression for — J2 a &Ai ^ y . ^ ^Q^( a ) (n^) + D^^j . 



3. Solution for the Electromagnetic Field 



With the current at hand, it is now possible to find the expression for the electromagnetic 
field as a function of the particles dynamical variables. However, note that in the expression 
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for j (x) the electromagnetic field appears everywhere through the Berry phases and explicitly 
through the magnetic field. This means that in fact reinserting j (x) in the expression for 
A (x) allows only to compute this last potential perturbatively. Assuming, as before that 
the current is relatively weak, the right hand side can be expanded at the first order in the 
field. The expressions depending on the field can be put to the left of the Maxwell equation, 
which can be now written : 

d 2 



[□ + 50 (x)] A (x,t) = / d 6 y 



5^5 6 (x - y) + 



1 



dx i dx^ An |x — y| 



|A=0 



with : 

% (x) 



■H, 



|A=0 



5 A (x) 



E 



dv (r^) 



dnf (x,t) 



~.R 



(si 



> A<^ 



+ A 



j,i\A=0 



)} H 

/J|A=0 



x — r 



Or 



(a) 



(«) 



OA 



dlin 



(a) 



,R 



(«), 



„(«) 



+ 4 



n,i\ A=0 



&4 



-iff 



(a) 



(r(«)) 



x — r 



J |A=0 



- £ [(„ (««■>) x V),i (x - ,<•>) + [A* (x"») x V 8l W (x - r<«>)]J ^ 

a 

The operator 50 (x) will in fact yields negligible contributions in our applications. We give 
nevertheless its form in the appendix for the sake of completeness. Note that the correction 
50 (x) to the Dalembertian are of order h 2 due to the definition of the terms involved. 
As a consequence, we can thus write for the electromagnetic field : 



A l U 



,t)= j G y (x 



with 



G^(x-y,t-t') 
1 



-VI- 50 (x) 



(x-z,t-t') 



dt 2 



g*S 3 (z - y) + 



d 2 



1 



dx l dxi Air |z — y| 



dz 



A more tractable for the effective Green function Gij (x — y, t — t') can be found if we move 
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to the Fourier transform. Actually, one can write : 



Gn (x - y, t - t') = J exp (ip. (x - y) - iu (t - t')) (s ir - ( cj 2 - (p 2 - 80 (p)) ) 



exp (ip. (x — y) — iu; (t — t')) 



co 2 -p 2 \ rj ' 2 CO 2 - p 2 



= y exp (ip. (x — y) — icu (t — t')) G-ij (p,u) dpduo 
and the operators 80, 8G are given in Fourier components by : 



X _ ViVj 



(p,w) 



<JOy ( P ) 



C<J 2 — p 2 

^> lrJ (p) 



G„(p^) = G«?(p^)(4, + i*^ 

\ Z W — p 

and the Fourier transform 80\ (p) is given in the appendix. 

Now that we have given the effective photon propagator, note that we can decompose 
the vector potential at order ti 2 at a point x : 

A (x, t) = J G tj (x - y, t - t') A(y, t')dy 3 dt' 

- [dt'J2 ( x(a) ) + ( x(Q) ) V fifc ] e lmi V m Gij (x - x< Q > (*') , i - i') dt 
= J Gij (x - y, t - t') ji(y, t')dyW - y dtMfe^VlGij (x - x<°> (0 , t - t') 
where .Mj is the effective magnetization differential operator : 



Mt = [Vm (X (Q) ) + ft (xW) V fifc ] 



A=0 



and the current ji is given by : 

9 (i a L + Dia) 



3i = eJ2 

a 

-E 



fcrf (x,t) 
^ ( r («)) 



- 2 R ' a) V r w + ^|A=o)] 5 (x - r™) 



dr- 
ee J 



(a) 



(^)- ^ - A 1 V p (a) + j4 n ,i|A=0 



&4 



(a) 



(r(«)) 



8 (x - r^) 



In recalls that the gradient is taken with respect to the variable x. 
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4- Electromagnetic part of the Hamiltonian 



We can now compute the Hamiltonian for the electromagnetic field as : 

dy 3 dx 3 A i (x,t) 



^-v 2 V — 

dt 2 I dx l dx'i 



A*(x,t) 



d 2 



dx l dxi 



G jm (x-z,t-t")r(z,t") 



dy 3 dx 3 dt'dt?(x, t)G (x - y,t - t') f (y , t') 



where the effective propagator G (x — y,t — t') is given through its Fourier transform : 
G (x — y,t — t')= / exp (ip. (x — y) — iuo (t — t')) G^ (p,iu) dpdu 



with : 



Gij (p,w) = (p,u) [(Si m (p 
150 ir (p) 



2 co 2 ) 



PlPm) ] G m j 



(p,cu) = <5 ir + 



5,: 



2p 2 -w 2 / p 2 — c<j 2 

1 50 ir ( P ) <^j 



(<W (p 2 - UJ 2 ) - pip m ) 



X _ PmP. 



p 2 — UJ 2 



s j -t 



2 p 2 - cj 2 



PrPj 



+ - 



1 <Ji. 



p 2 dG> s j (p) 



p 2 — a; 2 2 p 2 — w 2 p 2 — cj 2 2 p 2 — cj 2 p 2 — u 2 

If we neglect the retardation effect, which is legitimate since we work in the non relativistic 
limit, we can discard the dependence in w, so that : 

W 60 sj (P) 



Gij (P) 



+ 



+ 



and : 



^ / ^n 2 m + ^ | c? 3 a; (V x A 



W(x,t)G(x-y)j' (y,t) 



5. Electric potential 



Now, the case of the electrostatic potential is a bit peculiar. As we said before, we replace 
it through the potential equation in the Coulomb gauge : — V 2 A = J , so that 
/ d 3 xeA (x) Jo (x) = I J d 3 x J d 3 yj (x) i ^ J (y) = / £>x (VA ) 2 



57 



Considering v (r^) alone, one has first : 

v{r (a) )-^ J d 3 x(VA ) 2 
= J rf 3 xeA (x) e5 (x - r< a >) - l - J d 3 x (VA f 
= l - J d 3 xeA (x) e5 (x - r (a) ) 

^/ A / dV ° (x) ^b J ° (y) 



2 ^4vr rH-r^ 3 ) 



(a)_ r (0) 



) 



where V ^ r ( Q )— r (/3)^ i s the usual Coulomb interaction. The same token is applied to the 
derivatives of v (r^) appearing in the Hamiltonian, so that the part involving the electro- 
static potential can be written : 



h 2 



a 

[4°' W . A* 



EK 

.M(«) 



(a) 



Ar(«)ju"(«) 



V /m(«)jv(«) 



(a) 



jV"(«)M( a ) 



R 



A 



(«) 



V„(«)V 



J\f(a)jv(a) 



R (a) V R (a) 



( r (a)) 



£ 0M(°) £ 0AT(«) 



rf 3 x (VA ) 



4E 



+ E 



(«) 



E 

.M( Q ) V 

4 a) (x),a" 



iV(a)M( a ) 



^ 



)AT(«) 



NWM( 



A 



M(<*)AT(«) 



V r(q) V r ( 



We can now come back to our previous remark about the replacing the electrostatic po- 
tential before or after diagonalization. Had we chosen to replace this potential before ap- 
plying the diagonalization procedure would have in fact led to the same contribution to 
the Band Hamiltonian. Actually, it is straightforward to check that this early replacement 
would have led to consider the potential \ J2 a V (R^ — R^) in the Hamiltonian (and no 
more the potential v). This expression mixing the particles at an early stage would have 
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altered the diagonalization process. Actually, the diagonalization matrix wouldn't have 
been the simple tensor product of individual particles diagonalization matrices, but would 
have mixed the various spaces at first order in h. The result would be, at the order h 2 , in 
crossed contributions in the Berry phase for (not for R/ Q ) as can directly been checked). 
Namely, the Berry phase for would acquire a supplementary contribution proportional 
to V B (a) V p (/s) Es^a V ( R(a) - R(/3) ) ( A o™) S nWmW . Given our generalized results 

for the Diagonalized Hamiltonian, this term will not contribute to the Band Hamiltonian 
nor to the dynamical operators at order h 2 . This is the reason why we choose to treat the 
electrostatic potential in a less rigorous, but symmetric with respect to the vector potential, 
way. 

6. P particles Hamiltonian 
We can now gather all our results to obtain the final form for the Hamiltonian 



H 



N 



+ 



E 4% (*"") + \Y. V (* (a) - - \ E (*< (-""J • B<1 " + B< "' (* (<,) ) 1 

a. a. a. 

-fi k (x) .V Rfc B (Q) (x (Q) ) + B (Q) ./2.B (a) 



h 2 ^ 



EK 



.M( Q ) 



d(«) \ 

*/-\.n 



+ 



h 



T 



^ ( x ) ' -A) 



jV(«)M( a ) 

- Ar(«)M(«) 



v «!« V H!->E' / ( r(o, -'' m ) 



(a) 



V„(a)V 



p (a) V R (c) 



dyW^(x,t)G(x-y)f (y,t) 



D. Application : Hamiltonian at the lowest order 



We have seen in the previous section that, due to the complexity to the Maxwell equa- 
tion, the Full Hamiltonian to the order h 2 cannot be computed exactly but has in fact to be 
computed as a power series in the current. The Formula we gave included first order cor- 
rections for the Green function. Here, to give an application of our method we will neglect 
these first order corrections and consider the Hamiltonian at the lowest order in current, 
that is quadratic in current. We will see that it will result in a kind of generalized Breit 
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Hamiltonian, including magnetization-magnetization interaction. To do so, we divide our 
work in two steps. First computing the Hamiltonian without electromagnetic interaction 
and second including the electromagnetic corrections. 



1. Coulombian interaction 



If we neglect, the magnetic interactions, three simplifications arise. First, we can cancel 



R 



(«)/3 



all contributions to the magnetic field. Second, in the first order Berry phases, A l is 



R (<*) 

null, except if a — (3. As before, we will note A ' the zeroth order Berry phase, keeping 
in mind there are no crossed terms. Third, the crystal momentum dynamical variables are 
unchanged through the transformation. Thus = K^. The Hamiltonian reduces thus 
to : 



n 2 ^ 



h 2 ^ 



E K' 



A^ m ) 



M(a)N( a ) 



[4" ) M.-4o 



- jV( a )iVf( a ) 



R 



(a) 



o(«) V „(a) 



where the dynamical variables in that set up are given by : 



P («) _ 

1 N — 



» 



with 



AT( a ) 



+ 2 (W)iv(«)pw (• 4 o')p(«) JV (a) + (• 4 °')m««)p(««) (^pWatw) * 



2 

+H.C 



p(a) N (a) I °P(«)W( a ) - 



Vp^(R) 
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2. Electromagnetic interaction at the lowest order and magnetization-magnetization interaction 



We now turn to the electromagnetic part of the interaction in our Hamiltonian for P 
particles. It is found by isolating, at the second order in the current the interaction terms 
particles-electromagnetic field plus the internal Field Hamiltonian. 

In the non relativistic limit which is of interest for us, it takes the form : 



magnetic 



+ 2 6 



1 2 (9 2 (kW) + D&) + V (rW - rW>) 



|A=0 



|A=0 



- i ( xH ) - B(a) ( xW ) + BW ( xW ) ^ ( x )) " A* ( x ) - v ^ B(a) ( xW ) + B (a) ./2.B (Q) 

+ i y ^W^(x,t)G(x-y)f (y,t) 
Moreover, in first approximation, G (x — y) can be approximated by G (x — y) so that 
1 [ ^ 3 ^(x,t)G(x-y)f(y,t) ~ ± ^ dy 3 dx%(^ t)G# (x - y) £(y, t) 



rf 3 xd 3 yji(x,t) 



( x -y) l ( x -y) 7 



|x-y| 2 



2|x-y| 



As seen before, the current decomposes as : 

ii(y,t)=j li (y,t) + Mr / (y,t) e " Bi VW 

(we have used the fact that in the definition of the magnetization, the gradient of a delta 
acts as minus a derivative) where the magnetization M. e / and the current ji is given by : 

Mf (y) = J2 M *i f{a) ( x(a) ) 5 (y - r(a) ) 

= Eh( x(a) )+^( x(a) )v Rfc ]5(y-r (a) ) 

with .M^ <a ' ) is the particle a individual particle magnetization, and 

, , ^ 9 („<->) + fl<«>) 

jii(y,t) = e^^^- 



(a) 



fe-^ ) V r( . ) + ^)l 5(y-rW) 

V r j /J|A=0 



^ ( r («)) 



dr- 
at 3 



(a) 



•9^ 



ra,i|A=0^ 



+ 



dA i (r(«)) 



5 (y - r< a >) 



|A=0 
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ultimately, as seen in the previous section v (r^) can be replaced by J2p V ( r(o) - r(/3) )- 

These formula allow to derive directly the potential and the magnetic field in our approx- 
imation : 

2|x-rW| Ix-rWI 2 



b = x (x - rg) v>< ( r M) x (x - rjg 

^ |x-r( Q )| 3 Ix-rHI 3 

We can now compute the various terms appearing in -ff magn etic- The first term, involving the 
magnetization is given by : 

~ E (x (a) ) -B (x^) + B (x<«») ./i (x< a >)) - jl k (x^) .V Rfc B (x^) + B./2 (x^) .B (a) 

a 

= -\H (W ( r(Q) ) e^G?" (r<°>-r<«) ^(rM) + H.C. 

(rW) £^V fe V m G# (r<°>-r<0) jy(r^)) 
+ ]T £ Ip (r< a >)e' mi V m G# (rW-rt«)jy(rW) e P mi V fB Gj' (rW-rW) j^r^) 

a,P,a^P 

+ ^ ^(r( a ))^V mp V r V p G^'(r^-r^)/i m (r^) 

a,P,a^P 

The gradient, in all these expressions and in the sequel, is understood as acting on the first 
variable (here r^). 

The second relevant term is the energy of the internal magnetic field : 

l - J ^ 3 5 i (x ) t)Gnx-y)i i (y,t) 

= \ E M^gU^-^M^) 

+\ E Mr (a) ) + (r (a) ) 



a,P,a^P 
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We gather these two contributions to get : 

~\ E ( xW ) - B(Q) ( xW ) + BW ( xW ) ^ ( x )) " W - V ^ B(a) ( xW ) + B (a) ./2.B (Q) 



+ 1 1 ^W(x,t)G(x-y)f (y,t) 



\ E ^(r (Q) )^'(r (Q) -r (/3) )^(r (/3) ) 



+ ^i P (r (a) )e lmi V m G^ (r^-r^) 3lj (rWy™V m G% (r^-rW) ^(rW) 



The last term of interest for us, can be developer in the following way : 
'^tko ( kW ) + + (rW - rW) 



E< 



A j (r^) 



+ 2 6 



|A=0 



dkfdk^ 



A 3 ( r W) ^fc £ r (a)) 



|A=0 



i /^ 2 f4^)( k(a) ) + ^ (Q)N ) 

-^ Jl(r M)A(rW) + leM V ^ J tt) M (r' ">) .4>"") 



|A=0 



E Ji*(rW)G?(rW-r<«) Jni (rW) 



1 



2 
1 2 



£ (rW) + /if (rW) V*) £^V m G# (r<°>-r<0) ^ (r<«) + H.C. 



dkfdk[ a) 



a) 



|A=0 



The first equality is implied by construction of ^(r^) since this last quantity was precisely 
cleaned as E „e ^(^■..(^■)^'^^(^-^-) 



|A=0 
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We can now use the following facts : 

(x - y) = 
^ mi V m G^(x-y) = 

£ lri £ jm PVrVpG ij (x _ y) = 



S ij + 



( x ~y)i( x -y)j 

|x-y| a 



2|x-y| 

I |3 

l x -y| 

5 ij - 3 (x ~^ ( y "7 )j - f ^'5 (x - y) 



|x-y| 

Reintroducing the effective individual magnetization operator as : 

M ef(a) ( r (a)) = ^ ^(a)J + ^fe ^(a) j y fe 

one can rewrite ultimately at the order h 2 : 



if, 



magnetic 



= -E 



-(«)- r (/3) 



4e 



L 



e/(a) 



2 | r (a)_ r 09) 
( r W_ r 05)j 



| r («)_ r (/3) I 3 



xjiP) +H.C. 



4E 



M e f( a \M ef< >® - 3 



JH*).(r(«)-rM)jM < W).(r( l ')-rM) 



•(«)_ r (/3) 



_ ^L_A^e/(°)._A^e/(/3)^ ( r (a) _,.(/?)) 



| r (a)_ r (/3)| 



C r H_ r (/3)^ \ / f r (a) _ r 03)> 



f 



r («)_ r (/3) 



| r (a)_ r (/3)| 



+ E 





r (a)_ r (/3) 


2 
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Note here that the double scalar product involving the 2 tensor /^(r^) is performed on each 
of its indices separately. 
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The interpretation of the various terms can done by comparison to the very similar Breit 
Hamiltonian for the Dirac electron. Given the current is decomposed in two parts due to the 
velocity of the particles and their magnetization, the magnetic part of the Hamiltonian is 
mainly a current-current interaction. The first term is the usual current- current interaction 
involving ji(r^- a ') only, whereas the second term mixing ji{r^) and the magnetization of 
other particles is a magnetization-orbit coupling between different particles. This interaction 
is formally similar to the spin-orbit interaction term, except that here we are here in a non 
relativistic context, and that the magnetization arise as a band phenomenon. The third 
term is for the magnetization-magnetization interaction term which is formally similar to 
the spin-spin interaction term. The two last terms are new with respect to the usual Breit 
Hamiltonian. The fourth one is again a magnetization orbit coupling but of higher order, 
since it involves triples of particles. The last term is a correction to the energy due to the 
development to the second order of the free energy. It is a shift of the individual particles 
energy operator due to the field created by other particles. It should be present in the Breit 
Hamiltonian but is in fact neglected while considering the quantum field derivation of this 
last one since it is of order higher than for the Dirac case ( 301] ). 



3. Full Hamiltonian 



Since our development was in the second order in the currents we can simply gather 
the magnetic part and the electric part of the Hamiltonian to obtain ultimately the full P 
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particles Hamiltonian : 
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Recall at that point that we have performed a development in currents (or field) to the lowest 
order. As such in all our expression the electromagnetic field is set to 0. As a consequence, 
the variables k.( a \ are computed with A = 0. 
The dynamical variables are thus given by : 

-(«) _ r («)_|_ A nn 

L N — ^ Jf(.<*) 

i c ( Q ) _ k» 

JV 

Apart from the magnetic part of the Hamiltonian, already discussed, we have included the 
usual electrostatic part, as well as the Darwin term, which is analogous to the eponymous 
term in the Breit Hamiltonian. Note however that this term in our context is not only a 
contact term (i.e. of Dirac delta type) since its form depends ultimately on the form of the 
Berry phases that depend on the structure of the system at stake. 

Let us also remark ultimately that the expansion performed, is different from the one 
chosen in js(3] to derive the Breit Hamiltonian. Actually, the derivation of this last one 
through the context of quantum field involves a power expansion in — , which is in fact 
performed at the second order. 

V. CONCLUSION 

The diagonalization of the Hamiltonian for a Bloch electron in a magnetic field is an old 

n 

problem is solid state physics initiated principally by Blount |4j, who developed a general 
procedure for the removal of the interband matrix elements based on a asymptotic series 
expansion in the fields strength. In this way, the actual effective one-band Hamiltonian was 
obtained to the second order in the magnetic field. In this paper, we came back, in a way, to 
this old problem, but by trying instead to derive an in-band Hamiltonian as a series expansion 
in the Planck constant. A first attempt in this direction was done when we provided a 

n 

procedure at order H for an arbitrary matrix valued Hamiltonian [15]. This method resulted 
in an effective diagonal Hamiltonian in terms of gauge- covariant but noncanonical, actually 
noncommutative, coordinates. It has also revealed that a generalized Peierls substitution 
taking into account a Berry phase term must be considered for the semiclassical treatment of 



electrons in a magnetic field 25j . In particular, the Bohr-Sommerfeld quantization condition 



when reformulated in terms of the generalized Peierls substitution leads to a modification of 
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the semiclassical quantization rules as well as to a generalization of the cross-sectional area 
derived by Roth [5]. Soon later, the semiclassical diagonalization was extended to any order, 
when we developed another method of diagonalization where h was considered as a running 
parameter. This method allowed us to solve formally the problem in terms of a differential 
equation for the diagonal Hamiltonian with respect to Planck constant h, which could, in 
principle, be solved recursively by a series expansion in ft [l4]. As an example, the energy 
spectrum of a Bloch electron in an electric field was derived to all order in h. However, in 
the presence of a magnetic field, the differential approach beyond the semiclassical turned 
out to be very complicated. 

In this paper, we solved the Hamiltonian diagonalization for a Bloch electron in an 
electromagnetic field to the second order in h by applying a radically different procedure, 
which was recently developed and provides a general diagonal expression for any kind of 
matrix valued Hamiltonian 21[ (this diagonal expression is also solution of the differential 



equation of fl4|). The main advantage of our method with respect to Blount's procedure 
is obviously that our result is valid even in the presence of strong external electromagnetic 
fields. But a second important advantage with an expansion in h is that it allows us to 
consider particles in interaction. Indeed, although the removal of the interband matrix 
elements to the second order in h for a Bloch electron is important in itself, the principal 
objective of this paper was the determination of the effective Hamiltonian of interacting 
Bloch electrons living in different energy bands. As we have shown, even if the electrostatic 
interaction dominates the magnetic one, effects like magnetic moment-moment interaction 
mediated by the magnetic part of the full electromagnetic interaction requires a computation 
to second order in H. Comparison with other methods is difficult because to our knowledge 
other results for interacting Bloch electrons do not exist. Nevertheless because of the strong 
analogy with the Dirac equation it is possible to try a comparison with the Breit Hamiltonian 
for the Dirac electrons [23[. As for Breit, we found that the electronic current is made 
of two contributions: one comes from the velocity and the other is a magnetic moment 
current similar to the spin current for Dirac [22j ]. It is this last one which is responsible 
for the magnetic moment-moment interaction similar to the spin-spin interaction for Dirac. 
Another important interaction revealed by our approach is the moment-orbit coupling by 
analogy with the spin-orbit one. From the results of this paper, our goal in the future is to 
consider the physical relevance of interactions bringing into play magnetic moments. 
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VI. APPENDIX 



Since the right hand side of the Maxwell equation involves 
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where the operator 50\ (x) is obtained by isolating in — SA /^ y t ^ H rnat the linear term in the 
electromagnetic potential. Collecting the relevant terms yields directly : 
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The Fourier transform of this operator is straightforwardly obtained as being equal to : 
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